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Abstract 

We discuss various properties of the conformal field equations and their 
consequences for the asymptotic structure of space-times. 



o 

1 Introduction 

The causal structure, which decides whether a signal can be sent between two 
\ space-time points, is of fundamental importance for all of field theory. In the 

i«j '■ case of Einstein's theory of gravitation, general relativity, it acquires a particular 

^ | significance, because the null cone structure underlying the causal relations is 

itself a basic unknown of the field equations and as such subject to the causal 
relations in the small. 

From the beginning of general relativity the effect of 'deflection of light 
r> | in gravitational fields', demonstrated nowadays so clearly by the gravitational 

■ lensing phenomena, has been considered theoretically and observationally as one 

of the clearest manifestations of the deviation of the light cone field in general 
relativity from that of special relativity. As seen in the article by G. Galloway 
in this volume, which shows, among other things, the ubiquity of caustics on 
wave fronts in solutions to Einstein's field equations, it is not only a local effect 
but also related to the global structure of space-times. 

This relation to the global structure of space-time becomes quite pronounced 
in the case of black holes, which seem to be on the verge now of also becoming 
'observed' phenomena. The observational situation as well as various theoret- 
ical aspects of the idea are addressed directly in P. Chrusciel's article in this 
volume. However, tacitly or explicitly, the related questions of gravitational col- 
lapse, merger of black holes, generation of gravitational radiation, singularities 
of gravitational fields, etc. motivate to a large extent all studies in this volume. 

"To appear in: J. Frauendiener, H. Friedrich (eds.), "The Conformal Structure of Space- 
time: Geometry, Analysis, Numerics", Springer, Lecture Notes in Physics, Berlin, 2002. 
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The dual and critical role of the light cone structure in Einstein's theory 
and in physics in general led some researchers to speculate that the null cone 
structure might be a classical space-time reflection of a fundamental feature of 
a quantum theory of gravitation. This motivated efforts to recast the classical 
theory into forms which put the null cone structure, or associated structures 
like wave fronts, into the foreground. The article by S. Frittelli, N. Kamran, 
and E. Newman as well as the article by R. Penrose in this volume are aiming 
into this direction. Trying instead to derive qualitative and quantitative results 
about the solutions, the remaining articles explore and exploit, in one way or 
other, the light cone structure along more conventional lines. 

That in studies of the global structure of gravitational fields the light cone 
structure acquires a particular technical interest, has been demonstrated most 
impressively by the proofs of the singularity theorems ( f5l]| , |o7|| , cf. also |32|). It 
provides the most important tool for analyzing space-times in the large. While 
these theorems show that under reasonable assumptions space-times are singular 
in the sense that they contain non-complete causal geodesies, they do not provide 
any information on the specific nature of the singularities. It is expected that 
in general curvature singularities will arise, however, the theorems give not the 
slightest hint whether this will be the case, not to mention the causal nature of 
the singularities or their strength. 

In spite of various attempts to obtain more detailed results, the related 
problem of cosmic censorship (cf. 72 , still poses some of the most 



important outstanding questions of classical general relativity. While the limited 
use of the field equations made in the singularity theorems contributes to their 
generality, the same fact is responsible for the scarcity of details supplied by 
them. In view of the enormous richness of pathological behaviour observed on 
the level of pure space-time theory (cf. |2^] and the literature given there) it 
is not surprising that attacks at the problem of cosmic censorship which rely 
mainly on the topological methods used in the proof of the singularity theorems 
met little success. 

To disclose the more subtle secrets of Einstein's theory we need a deeper un- 
derstanding of the causal structure in the context of Einstein's field equations. 
In any serious investigation of the problem the consequences of these equations, 
which represent the principal restrictive principle of the theory, must be taken 
into account in all completeness. Getting control on the evolution process de- 
fined by Einstein's equations still remains the most important technical task in 
classical general relativity. 

There arises the question of how to combine the analysis of the field equa- 
tions with that of the causal structure. In Einstein's theory these structures are 
particularly closely related. By the local causality requirement the causal struc- 
ture is determined by the null cones of the metric. Conversely, provided certain 
pathologies are excluded, the causal structure can be used to reconstruct the 
null cone structure. Moreover, the physical characteristics of the field equations, 
which largely govern the evolution process, are given by the null hypersurfaces 
defined by the solutions. These tight relations between causal structure, null 
cone structure, and characteristics structure clearly offer possibilities for the 
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desired analysis, but it may be noted that they are also responsible for the 
quasi-linearity of the equations which, together with their gauge hyperbolicity, 
is the main source of our difficulties. 

The null cone structure is equivalent to the conformal structure, given by 
the information which is retained if the metric is subject to arbitrary conformal 
rescalings g — > Vt 2 g with conformal factors > 0. This may allow us to relate 
the null cone structure in a new way to the structure of the field equations and 
to the large scale behaviour of their solutions. 

However, to what extent can we draw from the behaviour of the field equa- 
tions under conformal rescalings of the metric conclusions on the large scale 
behaviour of the solution? The present article is trying to give some answer 
to this question. The latter might be rather aimless were it not for a concrete 
proposal pointing precisely into this direction. Motivated by earlier studies of 
the important physical question of how to characterize gravitational radiation 
without relying on approximations, R. Penrose suggested to analyses the asymp- 
totic behaviour of gravitational fields in terms of extensions of the conformal 
structure to null infinity (|3jJ, p6| ). 

The subsequent attempt ( fl23fl , pif) to resolve the controversy concerning 
the assumptions made in |i5| and to understand the smooth interplay between 
the conformal structure and the field equations suggested in |6(| disclosed the 
very peculiar behaviour of Einstein's field equations under conformal rescalings. 
They turned out to be conformally regular and assume with this property a very 
special position between conformally singular equations like the massive Klein- 
Gordon equations on a given background and conformally invariant equations 
like the Maxwell or Yang-Mills equations in four dimensions. It will be the 
purpose of this article to explain the origin, the properties, and the consequences 
of the resulting conformal field equations. 

Almost all the results described in the following have been obtained by con- 
sidering as a guiding problem the question to what extent the assumptions 
underlying ||| can be justified. We are beginning to see now that this question 
proved fruitful because these assumptions refer to the key geometrical struc- 
ture and draw such a sharp line between what is true and what isn't that they 
required, and still require, a most careful and detailed analysis. The results 
obtained so far indicate that in the most interesting case, that of a vanishing 
cosmological constant, the situation is much too subtle to simply allow us to 
say that the assumptions of |6j| are right or wrong. 

In the following we shall make a particular effort to sort out the delicate 
issues which arise here, hoping to put a controversy, whose 'arguments' too 
often relied more on personal taste than on insight, onto a rational basis. Such 
an effort is appropriate because due to the results which are available now 
the mathematical analysis is coming close to a point where the more subtle 
mathematical and physical questions concerning the foundations of the theory 
of isolated self gravitating systems become amenable to detailed discussions. 

The investigations are now way beyond a situation which would allow us 
to arrive at reliable statements by investigations of the linearized equations, 
'physical arguments' etc. In fact, it is difficult to develop any intuition on phys- 
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ical grounds for the extreme questions about solutions of Einstein's equations 
which are being asked here. Consequently, most results are stated in the form 
of existence theorems (referred to in the following only in a non-technical way). 
The physicist mainly interested in wave forms and observational facts, who con- 
ceives of this as merely interesting mathematics, may be reminded that abstract 
statements about the existence of general classes of solutions with certain global 
properties usually only epitomize collections of results about details of the solu- 
tions, particular features of the equations, and new techniques with many more 
applications. 

More importantly, it may be pointed out to him that the results alluded to 
above become of a practical use now, which was not anticipated in |]65| . As 
discussed in the articles by J. Frauendiener and S. Husa in this volume they 
form a basis for numerical calculations covering asymptotic parts of space-time. 
Most remarkably, though still to be explored in numerical work, they open 
the unique possibility to calculate numerically on finite grids maximal, globally 
hyperbolic, asymptotically flat solutions to Einstein's field equations, including 
their asymptotics in null directions and the radiation field at null infinity. These 
solutions will be determined solely and uniquely from standard Cauchy data 
(subject to certain, quite weak, conditions near space-like infinity) and require 
besides the idealization of asymptotic flatness no further approximation. 

Working out the details of such truly global solutions, either analytical or 
with numerical methods, will still require a considerable amount of work and 
inventiveness. However, the main point to be made here is that both possibilities 
are open now and both approaches have a lot to gain from each other. Each 
of them will pose questions which will shed new light on problems of the other 
approach and the answers obtained in one of them will suggest new solutions to 
the other one. 

A main open problem in this context is concerned with the question to what 
extent it will be possible to approach 'time-like infinity' within a numerical cal- 
culation if the solution develops singularities and event horizons. This problem 
has hardly been touched yet. Even for basic solutions like the Schwarzschild 
or the Kerr solution a careful and detailed analysis of time-like infinity in the 
conformal picture is not available yet. Some of these questions and also the 
related problem of gauge conditions will be addressed at the end of this article. 

It is clear that there is still a long way to go before we will be able to obtain 
wave forms of sufficient generality or arrive at insights which will bring us closer 
to a resolution of the problem of cosmic censorship. Nevertheless, it appears 
that the interaction between numerical and analytical methods referred to above 
offers the best prospects for further progress. 

Any work on general classes of solutions to the Einstein equations is nec- 
essarily of a very technical nature. In the following we shall give a survey on 
certain results and the current analytical work based on the conformal fields 
equations which emphasizes the general ideas. With the exception of the some- 
what detailed discussion of the conformal field equations, technicalities will be 
kept at a minimum. 
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2 The conformal field equations 



The core of Einstein's theory of gravitation is represented by his field equations 

Rfiu - A<^«, = k(T^ v - -Tg^), (1) 

for the metric , where A denotes the cosmological constant and T^ v the energy 
momentum tensor. Most properties of the gravitational field should follow from 
these equations. Those which don't and have to be 'put in by hand' characterize 
general features as well as specific properties of the physical situation to be 
modeled. For instance, in the case of cosmological space-times one may assume 
compact time slices while in the description of isolated gravitating systems one 
will require the existence of Cauchy hypersurfaces with asymptotically flat ends, 
but in both cases one has the freedom to prescribe initial data in agreement with 
these assumptions. In the following we will mainly be interested in the second 
of these scenarios. 

Equations (|l|) are known to be gauge hyperbolic in the sense that they are 
hyperbolic if suitable gauge conditions are imposed. Expressed in terms of 
coordinates satisfying the harmonicity condition V^V^x^ = 0, equations ([!]) 
take the form 

~g aP d a d p ~g^ = Q^(~g,d~g,T) (2) 

of a system of wave equations for the metric coefficients g M „, where the right 
hand side depends on the derivatives of these coefficients up to first order and 
on the energy momentum tensor. 

The form (|^) reveals the double role played by the light cone structure or, 
equivalently, the conformal structure. On the one hand the conformal struc- 
ture is part of the unknown since it determines the metric up to a positive 
point dependent factor, on the other hand it defines the characteristic cone 
{£ G T*M\ g af3 £ Q £/3 = 0} of the form (§) 01 Einstein's equations (for this and 
various other general notions and results concerning the field equations we refer 
the reader to (fH)). The characteristics strongly govern the local propagation 
properties of Einstein's equations. In fact, a more detailed argument involving 
implies that equations (Q) respect the local causality requirement of relativ- 
ity, provided the equations for the matter fields entering the energy momentum 
tensor are chosen appropriately. This is related to the fact that the characteristic 
cone determines the domains of uniqueness for equation ([l]) . 

The tight relations between the structure of characteristics, the causal struc- 
ture, and the light cone structure or, equivalently, the conformal structure, all of 
which depend on the solution at hand, is a very distinctive feature of Einstein's 
theory. However, the theory proposed by Einstein gives no reason to expect 
that the field equations would be particularly well behaved under conformal 
rescalings g^ v — > g^ v — fl 2 g^. The equations are chosen to determine isometry 
classes of solutions and are thus not conformally covariant. Formally, this is 
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seen from the transformation law of the Ricci tensor, given in n dimensions by 



R vp [g] = R vp [g] — V„ V p ft - g vp g xs - V A V 5 fi - — — V A n V 5 



where on the right hand side V denotes the connection associated with the 
metric g. 

If Einstein's equation (Q) is assumed in this relation and the conformal factor 
Q is considered as a given function on the solution manifold, we obtain an 
equation for the conformal metric g vp . It should be possible to analyse this 
equation by applying the techniques known for Einstein's equation, since the 
principal parts of the differential operators occurring in both equations, i.e. the 
parts containing the derivatives of highest order, are identical. For our intended 
application this equation for g vp has two deficiencies. The most conspicuous one 
is the occurrence of factors I/O, which makes the equation singular at the set 
J' = {il = 0}. Since this set will indicate some kind of infinity of the solution 
space-time, it will be of particular interest to us. The second problem arises 
because the function O, in fact also the manifold underlying the solution and 
the set J' , are not given a priori. They are related to the global geometry of 
the solution and must be determined jointly with the metric. 

Remarkably, Einstein's equations allow us to resolve these problems. They 
turn out to be conformally regular in the sense that there exist conformal rep- 
resentations of the equation which do not contain factors ft -1 (resp. factors fl 
in the principal part) and for which the conformal factor is determined by the 
equations themselves. To discuss this further we recall a few facts of conformal 
geometry. 

2.1 Conformal geometry 

A smooth conformal structure C on a manifold M is given by a family of smooth 
locally defined metrics whose domains of definition exhaust M and which are 
related on the intersections of these domains by conformal rescalings with pos- 
itive conformal factors. A smooth pseudo-Riemannian space (M, g) uniquely 
determines a smooth conformal structure on M which we denote by C g . Using a 
suitable partition of unity, it can be shown that any conformal structure can be 
obtained this way and we shall often represent a conformal structure in terms 
of a globally defined metric in the conformal class. 

A Lorentzian conformal structure C determines a null cone at each point of 
M. Conversely, the null cones fix C uniquely (cf. ]&| ). As a consequence 'space- 
and time-like' are well defined for C and a few further important space-time 
structures are implied. A null hypersurface for C is a smooth (say) hypersurface 
which at each of its points is tangent to the null cone field. The unique directions 
of tangency can be integrated to obtain null curves which rule the null hyper- 
surface. Considered as point sets, these curves coincide with null geodesies of 
the metrics in the conformal class and any null geodesic can locally be obtained 
this way. Extending the null geodesies generating a null hypersurface in any of 





(3) 
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their affine parameterizations we obtain an unique extension. The resulting set, 
called a wave front, will in general not any longer be a smooth hypersurface, 
because families of null geodesies tend to develop caustics, i.e. form envelopes 
or intersections. The wave fronts and their caustics are also invariants of the 
conformal structure. 

On the level of local differential geometry the most important invariant of a 
conformal structure C g is its conformal Weyl tensor C M vXp (observe the index 
positions), which is the trace- free part in the decomposition 

R" vXp = 2 {5 M [a L p]v - g v[X L p] ^} + C» vXp , (4) 

of the curvature tensor of g, where the trace parts are represented by the tensor 

V - ^ " 2(^1) R9 » v ) ■ (5) 

In terms of the Levi-Civita connection V of g and some coordinate basis d p our 
conventions are such that i? M v \ p X v = (Va V p — V p V \) X p . We denote by 
Ruv = R p upv and R the Ricci tensor and the Ricci scalar of g respectively. 

The class of torsion free covariant derivative operators naturally associated 
with a conformal class C comprises but is larger than the set of Levi-Civita 
operators for the metrics representing C. A Weyl connection for C g is a torsion 
free connection V on M which respects the conformal class C g in the following 
sense. For any C 1 curve x(X) defined on a neighbourhood of the origin of 
K and any frame field e^, k — l,...,n, which is parallelly transported along 
x(X) with respect to V, there exists a function S1a > along x(A) such that 
g(ei,ek)\ x (\) = £l\g{ei, ek)\ x (p)- R follows from this that 

V P # MI , = -2/pg^, (6) 

with some arbitrary 1-form / on M and that the connection coefficients defined 
by Va M d v = T M p v d p and Va„ d v = f M p v d p are related by 

T f , p u = T ll P v +S(f) IM P v , with S(f) fl P u = 6P ll f u +6P v f fl -g liU gP x h. (7) 

In the special case where the 1-form / is exact such that it can be written 
on suitable open sets U in the form 

f = -n- 1 dn, (8) 

with some smooth function Q > 0, the connection V coincides on U with the 
Levi-Civita connection V of the metric g satisfying g^ — fl 2 g pl/ . 

The curvature tensors of the connections V and V satisfying (Q) are related 

by 

& ux P - R» ux P = 2 { V [A S p] " v + S s " [A S p] 5 v ) = 2 {g» [p V X] f u 

+ V [p /" g X]v - g» v V [p / A] - g» [p f x] f v + g v[p f x] /" + g» [p g x]v fs f) , 
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were indices are raised or lowered with respect to g. This implies for the Ricci 
tensor R vp — W „ w and the Ricci scalar R — g vp R vp 

Ru P = Ru P - (n - 1) V p f v + V v f p + {n- 2) /„ f p - g vp (V A f x + (n - 2) fx f x ), 

R = R - 2 (n - 1) V A f x - (n - 1) (n - 2) fx / A . 

Setting here / = — f^ 1 cif2, we obtain for the metrics g, g which are related by 
the rescaling g = ft 2 g the relation (^) of their Ricci tensors, which implies for 
their Ricci scalars the relation 

Tt ' — 1 ~ n + 2 

4 — V ^ 6 - R[g]6 = -R\g]6^ . (9) 

n — 2 

For later convenience we set here = fl ~ . 

The analogue of the decomposition (Q) takes the form 

BP vXp = 2 {g» [A L p]u ~ g» „ L [Xp] - g v[x L p] p } + C vXp , (10) 

where the trace parts of the curvature tensor are represented by the tensor 

X f - n — 2 - 1 ,\ 

It is related to (||) by 

Vp f v — f p f v + - g^ fx f X — Lfj, v — L^. (12) 

By a straightforward calculation it can be verified that (||), ([n]) satisfy the 
identity 

Vp L pi/ — V^Z/pj/ = Vp Lp„ — Vp L pv + f x C jypp- (13) 
2.1.1 Conformal geodesies 

We have seen that null geodesies, considered as point sets, are invariants of a 
Lorentzian conformal structure. With any conformal structure Cg is associated 
in fact a much larger class of distinguished curves called conformal geodesies. 
We shall only give an outline of their properties here, referring the reader to 
[ p8| , p[ , and jl6| for details and further results. A conformal geodesic curve 
x(t) is obtained, together with a 1-form b — b(r) along the curve, as a solution 
to the system of equations 

(V ± xr + S(b)x%i X ± p = 0, (14) 

{^±b) v ~\b ll S{b)x ti v x x =Lxv± x , (15) 

where L and S are as in (||), (^). There are more conformal geodesies than 
geodesies for any metric in the conformal class because for given initial data 
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x* 6 M, x* € T Xt M, 6* € T* t M there exists a unique conformal geodesic x(t), 
6(t) near x* satisfying for given t* e R 

x(t*) = x*, :e(t*) = £*, &(t«) = 6*. (16) 

The sign of g(x, i) is preserved along a conformal geodesic but not its modulus. 
Conformal geodesies, considered as point sets, are in general different from the 
metric geodesies of a metric in the conformal class. Furthermore, they admit 
general fractional linear maps as parameter transformations. Initial data with 
g(x*,x*) — 0, x* ^ 0, determine null curves which coincide, as point sets, with 
the null geodesies of the metrics in the conformal class. The fact that we shall 
mainly make use of time-like conformal geodesies in the following, should not 
be taken as an indications that null or space-like conformal geodesies can not 
be useful. 

Conformal geodesies are conformal invariants in the following sense. Denote 
by / a smooth 1-form field. Then, if x(t), b(r) solve the conformal geodesies 
equations (|l4|), (|l5|), the pair x(t), b{r) — f\ x ( T ) solves the same equations with 
V replaced by the connection V = V + S(f) and L by L, i.e. x(t), in particular 
its parameter r, is independent of the Weyl connection in the conformal class 
which is used to write the equations. 

Assume that there is given a smooth congruence of conformal geodesies 
covering an open set U of our manifold such that the associated 1-forms b define 
a smooth field on U. Denote by V the torsion free connection on U which defines 



with the connection entering (14), ( |15| ) the difference tensor V — V = S(b) and 
denote by L the tensor ( |ll|) derived from V. Comparing with (|l2|), we find that 
equations (|l4]), (|l5| ) can be written 

Vi£ = 0, (17) 



L x „x x = Q. (18) 

Let efe be a frame field satisfying 

V ± e fc -0, (19) 

along the congruence, and suppose that 5* is a hypersurface transverse to the 
congruence which meets each of the curves exactly once and on which g(ei,ek) — 
©*?7ifc with some function 8* > 0. Following the argument which leads to (||) 
we see that g(ei, e^) = 2 rjik on U with a function satisfying 

V i e = 6<i,6>, 615 = 6*. (20) 

2.2 Derivation of the conformal field equations 

The curvature tensor always satisfies the Bianchi identity Vm i? M \y\\ P ] = 0. 
Cbserving (Q), we can write it in n > 3 dimensions in the form 

V [s C» HAp] = 2 {g v[x V S L p] »-g» [x V s L p]v }, (21) 
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which gives after a contraction 

V tl C»„ Xp = {n-?>){V x L pv -V p L Xu }. (22) 

We assume now that the physical metric g satisfies Einstein's field equation 
(Q) and try to express these equations in terms of the conformed metric 

g = Q 2 g, (23) 

and the conformal factor £1 > such that the difficulties pointed out above are 
avoided. Although our discussion can be generalized to include certain matter 
fields (cf. for instance |3^], |5j| and also, for an analysis with a different intention 
and different results, P. Tod, this volume), we assume for simplicity that matter 
fields are absent, i.e. T pv — 0. Equation (|22] ) reduces then to 

V p C" vXp = 0. (24) 
This equation is of particular interest for us because due to the identity 

v M (n 3 - n vXp ) = n 3 - n v M c v\pi 

which holds for any conformal factor, it has a certain conformal covariance. It 
implies for 

d» vXp = n 3 - n C» 1/Xp , (25) 

the Bianchi equation 

V M d" ux p = 0, (26) 

which, in a sense, represents the core of all the systems to be considered in the 
following. We consider now equation (^) for the conformal metric g and write 
it in the form 

R» uXp = 2 {g» [A L p]v - g v[x L p] + fi"- 3 d" vXp . (27) 

Using (^2|) also for g and observing (|25|) , (p6|), we obtain for L pv , which repre- 
sents the Ricci tensor of the conformal metric, the equation 

\7 X L p „-\7 p L x „ = n n - 4 \7 p nd^ Xp . (28) 



For later discussion we note how (|27|), fl28| ) generalize if we admit besides 
rescalings also transitions to Weyl connections which are related to the Levi- 
Civita connection V of g by difference tensors 

V-V = S(/), (29) 

with smooth 1-forms /. Since V - V = S^ 1 Vfi), it follows that V - V = 
S^ 1 d) with the smooth 1-form d p = Q f p + V^fl By (|Io|) equation @ is 
then replaced by 

R» vXp = 2 {g» [x L p]v - g» v L [Xp] - g v[x L p] + ft"" 3 d* vXp , (30) 
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while ( p8f ) and ( |l3| ) imply the equation 



V„L^-V |i i /W = n ^d^d^uxp. (31) 

By admitting rcscalings and transitions to Weyl connections we have ar- 
tificially introduced a conformal gauge freedom into the equations. If these 
equations are to be used to derive information about the physical metric g by 
analyzing boundary value problems, this freedom has to be removed again by 
imposing gauge conditions such that (i) the resulting equations ensure that <?, 
derived by ( |23| ) from fl and g, satisfies the field equations, (ii) interesting in- 
formation on g can be gained by suitably fixing CI and /. Two quite different 



methods have been discussed so far (|E3j, 24 , |35|) 



2.2.1 The metric conformal field equations 

We consider only conformal rescalings of the metric and try to get equations 
which allow us to specify the conformal factor in a useful way. If we observe (0) 
in the form R^ v — — R g pv in (||) , the resulting equation, which we read now as 
equation for $7, can be written in the form 

V IJ ,V v n = -QL fJ , v + sg fll/ , (32) 

with 

s = -V A V A »+ o 1 RSI. 
n L n (n — 1 J 

Equation ( |32"|) is obviously overdetermined. If we apply V p to both sides, com- 
mute V p with Vp, contract over p and is, and use the twice contracted Bianchi 
identity, we obtain the integrability condition 

V^s = -V / fii 1// x, (33) 

which we read as an equation for s. 

To see the role of the cosmological constant, we observe that the transfor- 
mation law of the Ricci scalar gives 

A = (n- 1) (2Qs- V p riV p n). (34) 

Equations ([32|), (|3^) then imply for the quantity on the right hand side 

- v A (2 n s - Vp n v p n) = v A n s + n v A s - v n v p v A o 
= v A n s + n (-v 17 n x, yA ) - v p o (-n l pA + 9pX a ) = o, 

which reflects the well known fact that the twice contracted Bianchi identity 
implies that A must be constant. It follows that ([m]) will be a consequence of 
the other equations if it is arranged to hold at one point by a suitable choice of 
the initial data. 
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We consider (§|), @, @, ©, (f|), (§4|) as equations for the unknown 
tensor fields 

and refer to them as the metric conformal field equations. These equations hold 
for any vacuum solution g and any conformal factor f2. 

A way to impose a gauge condition on the conformal factor is suggested by 
the following observation. If g is any metric and R* a given function, the metric 
g* = O 7 ^ g will have by (^) the Ricci scalar R[g*] = R*, if solves the equation 

7? — 1 n + 2 

4 — - - V M W 9 - R[g] 9 = -R* 9~ . 

Since locally this equation can always be solved for a positive function 9, it 
follows that locally the Ricci scalar R of the conformal metric can be prescribed 
in the metric conformal field equations as a gauge source function for the con- 
formal scaling. Its choice will depend on the given situation and the choices of 
the other gauge dependent quantities. We can set R = but, depending on the 
circumstances, there may be better choices. In any case this leaves the freedom 
to specify suitable boundary data for the conformal factor. 

From the way we obtained the metric conformal field equations it is clear 
that they are equivalent to the vacuum field equations (ignoring questions of 
smoothness) where f2 is positive but that they make sense also where fl van- 
ishes. While the system is regular in this sense, it is also complicated and highly 
overdetermined. We shall see below that the formal regularity of the equations 
does not necessarily imply that the data must be regular nor does it ensure by 
itself that the asymptotic behaviour of the solutions can be controlled conve- 
niently. To get an insight into this one needs to work out suitable initial value 
problems in great detail and control the constraints as well as the evolution. We 
now take a look at the various parts of the system. 

2.2.2 The Bianchi equation 

For a solution g to the vacuum field equations the Bianchi identity (^lj) takes 
the form 

V 5 C vXp + V p C v s\ + V A C» vp5 = 0. (35) 
Contracting with S7s, commuting derivatives, and observing (|24|), we obtain 

v« v* C" VTp = 2 {cr 

Thus, in the Lorentzian case, which will be considered in the following, we find 
the well known fact that the Bianchi identity implies a system of hyperbolic 
equations for the conformal Weyl tensor in any dimension n > 4. 

In n = 4 dimensions we can define the left and right dual of the confor- 
mal Weyl tensor, given by *C plJTp = \ e M „ 5n C S7TTp and C* urp = \ i Tp 5n C pv ^ 
respectively, and it is well known that they coincide. It follows that 

e-K ' p Vj „ r p = 2 V<5 C ^ v -n = 2 V5 W7r " = e^ 1 v p V5 C nT pj 
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which implies that the contracted Bianchi identity (EJ) is equivalent to the 



Bianchi identity (35). It follows that ( |24| ) and, by the same argument equation 
(p6|), imply hyperbolic equations. 

In dimensions n > 5 the situation is quite different as can be seen by a simple 
function counting. The conformal Weyl tensor has in n dimensions 

-^n 2 (n 2 -l)-^n(n + l), 

independent components while the contracted Bianchi identity ( pi| ) provides 

1 



— n in 



2) (n + 2), 



independent equations. Moreover, if we express (|2J) in Gauss coordinates with 
x° the affine parameter on the geodesies orthogonal to the space-like hypersur- 
faces {x° = const.}, it follows that the equations V«5 C s ox P — involve only 
derivatives in the directions d x i, . . ., d x n-i and thus provide 

1 . 

-n(n-l) 

constraints. Consequently, there are not enough equations in ( p^ ) if n > 6 
while for n = 5 there are as many equations as unknowns but not sufficiently 
many evolution equations. It follows that a system of conformal equations which 
includes as equation for the rescaled conformal Weyl tensor only the Bianchi 
equation provides a good evolution system only in four dimensions. 

If we use in n > 5 dimensions the full Bianchi identities to derive a sys- 
tem of wave equations for d 1 ^ v \ pi it appears that not all terms of the form 
f2 _1 V M n can be removed from the equations. It would be interesting to know 
whether there exist other conformally covariant equations which could replace 
the Bianchi equation in higher dimensions. If it turns out that a conformally 
regular representation for the Einstein equation can only be obtained in four 
dimensions, it would also be interesting to know whether there is a principal dif- 
ference in the asymptotic behaviour of the fields in four and higher dimensional 
space-times. 

We shall restrict our discussion in the following to the case of four dimen- 
sions. 

Instead of expressing the intrinsic hyperbolicity of the Bianchi equation in 
terms of a system of wave equations we work directly with the first order system. 
The benefits of this will be seen later. There are known two quite different 
possibilities to make the hyperbolicity of the equations manifest, both of them 
have been used in numerical calculations (cf. J. Frauendiener and S. Husa, this 
volume) . It has been shown in |]37f how to obtain hyperbolic equations from the 
Bianchi equation in the ADM representation of the metric. The other method, 
which is based on a frame formalism, has many advantages and can be adapted 
to various geometrical situations. We can use real orthonormal frames (cf. KJ) 
but most convenient are pseudo-orthonormal frames or spin frames. In the spin 
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frame formalism the equations split almost automatically into constraints and 



(symmetric hyperbolic) evolution equations (cf. p|, ]27| ], 1 32 ) . We shall use 
this fact here to discuss certain features of the Bianchi equation. 

In the spin frame formalism we have a frame field e aa ', o-,a' = 0,1, with 
eoo': en' real null vector fields and eoi', eio' complex vector fields which are 
conjugate to each other. Their only non-vanishing scalar products are given by 
g(eoo', en') = 1, g(eoi',eio<) = —1. In the following discussion it is useful to 
think of these vector fields as being obtained in the form 

1 / ^ 1 / 

eoo' = (eo + e 3 ), ew = —j= (eo - e 3 ), 

1 / • > 1 / • ^ 

eoi' = — ^ (ei + t e 2 ), e w > = —j= [e\ - i e 2 ), 

from a frame field e&, k = 0, . . . , 3, satisfying g{ei, e^) = r)ik- 

The rescaled conformal Weyl tensor is represented by a completely symmet- 
ric spinor field <fiabcd, and the Bianchi equation by the equation 

= A abca > = V f a - fiabcf, (36) 

where V aa ' denotes the covariant derivative on spinor fields in the direction 
of e aa ' and the index is raised with the antisymmetric spinor field e ab which 
satisfies e 01 = 1. 

If we set A Qbcd = A Qbca / r d a ' with r a a ' = S a 1 6 > a ' - S a 8 V °' we get the 
decomposition A abcc ; = A( abcd ) — | e^fc A a 6)/ * into its irreducible parts 

~2 A( abc(J ) = P (pabcd ~ 2 P( d * </> abc )/ , A 6/ * = 2? 6 ^ </> abe / ■ 

Here P — r aa V aa ' = y/2 V eo and 2? ab = T( a V b ) a / denote covariant directional 
derivatives such that 2? o = — Voi', An = 2?io = ^ V e3 , = ^10' ( CI - lUli 
for more details of the underlying space-spinor formalism). 

In a Cauchy problem one will in general assume eo to be the future directed 
normal to the initial hypersurface S. The operators D ab then involve only 
differentiation in directions tangent to S and the equations A ab y^ = are 
interior equations on S. They represent the six real constraint equations implied 
on S by the Bianchi equation. 

The remaining equations for <j) a bcd should be evolution equations. Multiply- 
ing the equations by suitable binomial coefficients (and considering the frame 
field and the connection coefficients as given), we find that the system 



4 

b + 



A( Qbc d) = 0, (37) 



is in fact symmetric hyperbolic, i.e. it has the form 

A^d^u = H(x,u), 
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with a C^-valued unknown u, a C w -valued function H (x, u) , and N x iV-matrix- 
valued, possibly u-dependent functions A 11 which are hermitian, i.e. T = A^ , 
and for which there exists at each point x a covector £ M such that A^^ is positive 
definite. There exists a well developed theory for such systems (cf. p8[ , j58j, 
[ [77| , and also H.-O. Kreiss and O. Ortiz, this volume). 

While the constraints implied on a given space-like hypersurface are deter- 
mined uniquely, there is a large freedom to select useful evolution systems. In 
the present case it turns out that any system of the form 

= 2 a A(oooi'): 
= (c - d) A( 00 ii/) - 2 a A (0000 , ) , 
= (c+ d) A (0U1 ,) - (c - d) A( 00 io')i 
= 2eA (llir) - (c + d) A(ono'), 
= -2eA (m( y), 

with a, c, e > and — (2e + c) < d < 2a + c, is symmetric hyperbolic (the 
system (j3?]) occurs here as a special case). This flexibility allows the Bianchi 
equation to be adapted to various geometrical situations such as characteristic 
initial value problems (Q, @, |§, ©)' standard Cauchy problems 

(Q, H|), and initial boundary value problems (|3§], @). Note that for 
these systems the set of characteristics, given by the hypersurfaces of the form 
{/ = 0, df ^ 0} with / G C^Af) on which det{A^ f^) = is satisfied at each 
of their points, does not only comprise the null hypersurfaces and can in fact 
be quite complicated. 

2.2.3 The constraints 

Related to the gauge freedom is the fact that, written in suitably adapted coor- 
dinates (and frame field, if used), Einstein's equations as well as the metric con- 
formal field equations imply inner equations, i.e. constraints, on any time-like, 
null, or space-like hypersurface. We will consider in the following the latter two 
cases, because they are particularly important in initial value problems. (Time- 
like hypersurfaces are of interest in initial-boundary value problems, however, 
even in general analytical discussions of this problem a detailed discussion of 
the constraints on these hypersurfaces is not necessarily required H].) 

Null hypersurfaces representing the physical characteristics of the Einstein 
equations, the system of inner equations is quite different in nature from those 
arising in the non-null case. There are on the one hand the intrinsic equations 
arising because of the gauge freedom and on the other hand the characteristic 
inner equations arising due to the intrinsic wave equation character of the field 
equations and these equations combine in a very peculiar way. Let us consider 
a situation were data are to be prescribed on a null hypersurface N and another 
hypersurface which has a space-like intersection E with N. Then some of the 
unknowns can be prescribed completely freely on N. All the remaining unknowns 
can then be calculated on iV in suitably adapted coordinates from certain data 
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on E by solving a certain hierarchy of inner equations on N, each member of 
which reduces to a system of ODEs along the null generators at their subsequent 
integration. Those equations which are not used in this procedure will then be 
satisfied on N if they hold on S (cf. Q for the vacuum case). 

The enormous simplification of the equations resulting from these facts in 
a suitably adapted gauge and the possibility to integrate the equations along 
characteristics 'out to infinity' led a number of groups to develop very efficient 
characteristic numerical computer codes (cf. R. Bartnik and A. Norton, J. Font, 
L. Lehner, this volume). 

Though being more complicated, the inner equations induced on null hyper- 
surfaces by the metric conformal field equations have a similar structure and, 
besides the data required to fix the conformal gauge, the free data are, of course, 
up to a rescaling the same ( Q ) . 

In the following we will mainly be concerned with Cauchy problems and 
therefore consider now space-like initial hypersurfaces. The standard constraints 
induced on space- like hypersurfaces by the vacuum field equations form a system 
of underdetermined elliptic quasi-linear equations for the data induced by the 
4-metric on the initial hypersurface, i.e. the inner metric h a p and the extrinsic 
curvature Xa/3- There exists a detailed theory for these equations which, under 
the assumption of constant mean extrinsic curvature, reduces the problem of 
solving the constraints to a problem of solving a linear elliptic system and a 
semi-linear elliptic scalar equation (cf . Q , [jl3| ) . 

The metric conformal field equations imply a system of conformal constraints 
(pt|, p5fl) on a space-like hypersurface S which comprises the usual vacuum 
constraints but which is more complicated. This is not only due to the artificial 
scaling freedom which has been introduced into the vacuum equations but also 
because the equations include integrability conditions and the Gauss-Codazzi 
equations. There appear to be two ways of dealing with these equations. The 
first method solves essentially the standard vacuum constraints and calculates 
from the vacuum data the conformal data after introducing a suitable conformal 
rescaling. This method, which has been analyzed with a certain completeness 
and underlies the present numerical calculations, introduces certain numerical 
difficulties (cf. J. Frauendiener and S. Husa, this volume). A second possibility 
would be to solve the conformal constraints directly. First steps in a systematic 
investigation of this approach are reported in the article by A. Butscher in this 
volume (cf. also the discussion in |j33f ). 

Any of these approaches has to deal with a problem which is not considered 
in the usual discussions of the standard vacuum constraints. In the situations 
which are most interesting for our applications the initial hypersurface contains 
a subset E, not belonging to the physical space-time, where the conformal factor 
vanishes. As seen in detail in the next sections, this set can be a 2-dimensional 
boundary which represents points at null infinity (cf. L. Andersson, this volume) 
or consist of a number of discrete compactification points, E = {ii,«2; ■ ■ -iiv}; 
representing space-like infinities (cf. S. Dain, this volume). In principle, E 
could consist of a union of such sets or of more complicated sets. In any of these 
situations there arises the question about the precise smoothness of the data 
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near E in suitably chosen coordinates. 

It turns out that in general, even under the strongest smoothness assump- 
tions on the 'free' data, there may occur terms of the form r k log' r near E, 
where r denotes the (non-physical) distance from E. However, it also turns out 
that under mild additional assumptions logarithmic terms will be absent near 
E. We will discuss some of the subtle questions arising here later. 

2.2.4 Gauge conditions, hyperbolicity, and local evolution 

In principle everything known about Einstein's equations also applies to our sys- 
tem. If we introduce a fictitious energy-momentum tensor by writing nT pv — 
2 {La V — Lg^y) with a constant k^O (introduced here to emphasize the simi- 
larity, we could set k = 1), our equations imply after a contraction of ( |2^ ) the 
system 

^[\Tp] v - g V[ A T g p ] u = - V^fitf v \ p , 
V M d" v x P = 0, 

K 1 

V M V„f2 = -- (T M „ - - Tg ttv )Q, + sg^, 

V l *s = -'^(T^-^Tg^)V v n, 

which can be read as Einstein equations with certain matter fields, in which 
R = —kT represents a gauge source function. We note that in boundary value 
problems for these equations it has to be made sure, by the way the problems 
are posed and the data are given, that equation (|27]), which gives geometrical 
meaning to the fields L pu and v \ p , will be satisfied. 

Following the discussion of the Bianchi equation above (or that of |37]], if 
one prefers to consider the metric instead of a frame field as a basic variable) , it 
is not difficult to extract symmetric hyperbolic evolution equations, the reduced 
equations, from the last four subsystems. Of course, any choice here will leave us 
with constraint equations. It remains to discuss the first equation. In principle 
one can bring to bear here all the methods developed for the ADM equations 
if one feels most comfortable with these equations. However, it should be kept 
in mind that this is only one of the many possibilities of dealing with equation 
(p7|). In the references given in this article there can be found many others 
which have been employed successfully under various geometrical assumptions. 

If one is dealing with the Einstein equation above or directly with equation 
( p7| ) , in either case one will have to make the decision whether to use the metric 
coefficients g^ or frame coefficients k as basic variables. Furthermore, one 
needs to choose gauge conditions, i.e. a way to specify the coordinates and, 
in the second case, the frame field. Usually this is done in terms of a gauge 
system of differential equations and suitable initial conditions. A large set of 
possibilities is known by now. 
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Basic objectives in choosing the gauge are the following (cf. also Q): (i) 
The gauge should be adapted to the type of problem (Cauchy problem, initial- 
boundary value problem, characteristic initial value problem), (ii) The evolution 
equations should be simple, manifestly hyperbolic, or of some other form (e.g. 
'elliptic-hyperbolic') which allows us to handle them in existence proofs or in 
numerical calculations (cf. the discussion by H.-O. Kreiss and O. Ortiz, this 
volume). There are available gauge condition which are governed by implicit 
wave equations like the harmonic gauge or the conditions for frame or spin 
frame fields employed in j27j, J|^], and there are known explicit conditions, in 
particular certain geometric gauge conditions, which lead to a reduction of the 
number of the unknowns (cf. Q , |37j and the discussion below) and which allow 
us to control to some extent the development of a degeneracy of the gauge, (iii) 
The gauge conditions as well as the constraints should be preserved during the 
evolution, (iv) The gauge should allow us to control its behaviour during the 
evolution in numerical or in analytical studies, (v) It should have a sufficiently 
long 'life time'. 

While (ii) and (iii) are in the case of hyperbolic equations local properties 
which are accessible to direct calculations, (iv) and (v) depend very much on the 
given situation and cannot be handled in a general way. These points belong to 
the most difficult ones in the large scale evolution problem in numerical as well 
as in analytical studies. 

There are known now various ways to perform hyperbolic reductions (cf. p7|) 
for initial value problems for the metric conformal field equations: given a suit- 
able solution to the constraints, we find a reduced system of hyperbolic evolution 
equations for which the local existence and uniqueness of solutions follow from 
standard theorems, and, by invoking the available integrability conditions, wc 
can derive a subsidiary system which allows us to show that all constraints and 
gauge conditions are preserved during the evolution if they are satisfied initially. 
The significance of the conformal field equations is determined to a large extent 
by the fact that they admit hyperbolic reductions irrespective of the sign of the 
conformal factor. 

After we use the large freedom to choose the gauge system and the reduced 
system, the subsidiary system is essentially fixed. For the analytical discussions 
of solutions, the latter is only required to possess the uniqueness property, which 
ensures that trivial data imply the trivial solution. However, since in the numer- 
ical treatment of the Einstein equations the constraints are only satisfied up to 
a numerical error, the detailed propagation properties of the subsidiary system 
may become important in that context. A systematic analysis of how changes 
in the gauge system and the reduced system affect the propagation properties 
of the subsidiary system and whether this dependence can be used to choose 
reduced systems with improved evolution of the constraints has not been given 
yet. 
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2.2.5 The general conformal field equations 

We will now consider a system of conformal field equations which employs the 
additional freedom gained by admitting general Weyl connections. Combined 
with suitable gauge conditions it will lead to a system of evolution equations 
which is quite different from and in fact much simpler than the evolution systems 
obtained from the metric conformal field equations. 

Assume that g solves Einstein's vacuum field equations with cosmological 
constant A and let g = O 2 g be given a conformal metric where O is some 
conformal factor. Furthermore, let b be a smooth 1-form. We denote by V 
and V the Levi-Civita connection of g and g respectively and by V the Weyl 
connection for g which satisfies V — V = S(b). It follows that V = V + S(f) 
with / = b — 9 _1 V0. We define a further 1-form by setting = 6 / M + V M 6. 

Let efc be a frame field satisfying g(ei, ek) — r\ik- We denote by V/c and 
the covariant derivative in the direction of ek with respect to V and V and define 
the connection coefficients 3 k of V in this frame field by =T i 3 k^k- If 
we express now all tensor fields (except the ek) as well as equations (|30|), (pl[), 
( p6[ ) and the condition that V be torsion free in terms of the frame field and 
the connection coefficients, we get for the unknown 

u={e^ k , f^' fe , L jk , (Tjki), (38) 

the system of equations 

[e p , e q ] — (Tp q — T q p ) e;, (39) 



e p(r<? j) e q{^p j) Ffc j(Tp q fg p) + Tp fcfg j Tq fefp j (40) 

= 2 {g l [ p L q]j - g 1 3 L [pq] - g ][p L q] l } + Qd r jpq , 

-bqj ^ ' q ^pj — did jpq-, (^l) 

Vid? jk i = 0. (42) 

In the last equation the connection V is used. This poses no problem because 
the connection coefficients of V in the frame ek are related to those of V by 
the purely algebraic relation IV k = r, k + S j j fk + 5 j k fi - Vik rf l fi, whence 
f t = l/nt t k k . 

No differential equations are given above for the fields 8 and dk = Q fk + 
VfcO, they reflect the artificially introduced conformal gauge freedom. Note 
that they cannot be prescribed quite arbitrarily. For a solution which extends 
smoothly into a domain containing a set where — 0, we should have d = V0 
on this set. 



To obtain the metric conformal field equations, (26), (|27|), (|28|) have been 
complemented by equations for the conformal factor which contain basic in- 
formation on Einstein's equations. Finding a suitable gauge to obtain useful 
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evolution equations was then considered a different issue. In the case of the 
equations above we will combine these two steps into one. The extended gauge 
freedom will be used here in a very special way. The question whether there are 
other useful possibilities has not been investigated. 

We begin by describing the construction of a gauge. Let S be a space-like 
hypersurface in the given vacuum solution (M, g). We choose on S a positive 
'conformal factor' 0*, a frame field ek*, and a 1-form b* such that g(e^, ek*) = 
0* Vik and eo* is orthogonal to 5. Then there exists through each point x* 6 S 
a unique conformal geodesic x(t), b(r) with r = on S which satisfies there 
the initial conditions x = eo* , 6 = 6*. These curves define a smooth caustic free 
congruence in a neighbourhood U of S if all data are smooth. Furthermore, b 
defines a smooth 1-form on U which supplies a Weyl connection V as described 
above. A smooth frame field ek and conformal factor O are then obtained on 
U by solving ([n]), ( p0| ) for given initial conditions ek — ek*, b = b* on S. 
The frame field is orthonormal for the metric g = 2 <?, i.e. gik = g(ei,ek) — 
diag(l, —1,-1,-1). Dragging along local coordinates x a , a — 1, 2,3 on S with 
the congruence and setting x° = r we obtain a coordinate system. In this gauge 
we have on U 

x = e = 8 T , V* = 0, L ok = 0. (43) 

Coordinates, a frame field, and a conformal gauge as above will be referred to 
as a conformal Gauss system. 

We could now derive from (|2(]) and the conformal geodesies equations an 
equation for dk = bk — O fk + VfeO and use it together with ( |2p| ) to comple- 
ment the system (|39| ) to (0). Surprisingly, however, the fields 9 and dk can be 
determined explicitly ( |35|| ): If g is a solution to Einstein's vacuum equations 
with cosmological constant \, the fields O and dk are given in our gauge by the 
explicit expressions 



T 



2 



A 1 



e = e* (l + r <6*,i*>+ y \ Q-^- + -g\b„K)) ) . (44) 

d = 9, d a =< 6*,9*e a * >, a =1,2, 3, (45) 

where the quantities with a subscript star are considered as constant along the 
conformal geodesies and 5" denotes the contravariant version of g. 



With these expressions (py), (40), (41), (42) provide a complete system for 



u. In spite of the fact that we use a special gauge, we refer to these equations 
as the general conformal field equations to indicate that they employ the full 
gauge freedom preserving a conformal structure. 

One might wonder why in the present case it should be sufficient to give the 
explicit expression for and 6 while in the case of the metric conformal field 
equation we had to use Einstein's equations again to obtain equations for the 
conformal factor. However, the occurrence of A shows that information about 
Einstein's equations is encoded in (Q), (f45|). We note the following important 
consequence. If we have a sufficiently smooth solution to the metric conformal 
field equations in a domain comprising a non-empty set = {Q = 0, d£l Q}, 
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it follows from ( |34| ) that V p Q V p il = —A/3 on J", i.e. the causal nature of J 
is determined by the field equations (|]65|). The analogous result 

V fc 9V fc e = -iA where 6 = 0, (46) 

follows in the present case from (|44| ) , ( flBf ) and the definition of b if the congruence 
does not degenerate where 6 vanishes (|35|l). 

Our gauge is not only distinguished by the fact that it is provided by the 
conformal structure itself and the explicit information on 6 and b, but also by 
a remarkable simplicity of the resulting evolution equations. Setting p = and 
observing the gauge conditions (^) we obtain from ([59]) - (0) 

d T e% = -f,' oe^i, 

9 T r g ! j — — Tfc 1 j r 9 k + 9 % ^gj + ff* j — SjO £ 9 * + 6 cf jOq , 

f) f A A% 

^qj u jOqi 

VjeP jki = 0. 

Extracting from the Bianchi equation by the methods discussed above a sym- 
metric hyperbolic system, we get symmetric hyperbolic reduced equations for 
those components of u which are not determined explicitly by the gauge condi- 
tions. It has been shown in |}5| for such a choice of reduced equations (assuming 
(p4|), (jiq)) that any solution which satisfies the constraints on a suitable space- 
like hypersurface does indeed satisfy the complete set of field equations in the 
part of the domain of dependence of the initial data set where 6 is positive. 

A system of similar simplicity is obtained if we use metric Gauss coordinates 
and the vacuum Bianchi identity to solve the vacuum field equations directly 
(p^j). While it is well known, however, that in general Gauss coordinates 
tend to develop caustics quickly, it will be seen below that conformal Gauss 
coordinates can be arranged to cover large space-time domains or even entire 
maximal globally hyperbolic space-times in non-trivial cases. 

Our main purpose in reformulating the field equations was to obtain re- 
sults about their solutions. In fact, all the properties above have been used to 
establish properties of solutions by formulating and analyzing boundary value 
problems which take into account as far as possible the geometric nature of the 
situations under study. In the following we shall indicate the kind of questions 
which have been asked and the type of results obtained so far. 

3 The Penrose proposal 

Experience with other theories governed by hyperbolic field equations suggests 
that Einstein's theory should admit wave phenomena. For the same reason for 
which some work is required to bring to light the hyperbolic nature of Einstein's 
equation, we need to make an effort to exhibit the wave phenomena. Since the 
theory only supplies the metric and the associated curvature field, we need to 
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explain which part of the gravitational field shows wave-like behaviour with re- 
spect to which reference structure and to what extent wave phenomena observed 
in mathematical expressions are of any physical significance. The standard text- 
book discussion of gravitational radiation bypasses this problem by postulating 
some ad hoc background solution and studying solutions to the linearized equa- 
tions on it. This method allows one to apply the machinery developed in other 
areas of field theory, but specific features of Einstein's theory arising from its 
non-linearity are suppressed. Nowadays there exist much more sophisticated ap- 
proximation methods, which also take into account insights from the approach 
to be described below (cf. @, f§). 

Work initiated in the late 1950's by Pirani, Trautman, Bondi and others led 
to a concept of radiation which involves an 'approximation to nature', namely 
the idealization of an isolated system, but which does not require for its defi- 
nition any mathematical approximations in the given framework. Clearly, the 
gravitational field can be expected to develop a wave-like character only suf- 
ficiently far away from the generating system. The idea then is to consider a 
gravitating system, consisting of stars, black holes, or bounded systems of stars 
like galaxies, which is so far away from other such systems that the latter can 
essentially be ignored and taken into account only in terms of the gravitational 
radiation emitted by them. Thus one is led to study the far field of this system. 

Following a light ray which escapes from the system to infinity we may 
expect the field to become weaker and weaker and in fact asymptotically flat 
in the sense that the metric field approaches that of a limit Minkowski space. 
If some substructure, e.g. some component of the Weyl curvature, implies a 
non-trivial field on this limit space it might attain a natural interpretation as 
the radiation escaping from the system. 

To convert this idea into a workable definition two things are needed: a 
precise description of the limit procedure and detailed information on the decay 
of the metric and the curvature field. Bondi et al. || , Newman and Penrose |3lJ , 
and Sachs |7j] specify the limits in terms of so-called Bondi coordinates, which 
are generated by distinguished families of outgoing null hypersurfaces which 
extend in the future smoothly to infinity. On the hypersurfaces a convenient 
parameter is specified on the generating null geodesies which themselves are 
parameterized by spherical coordinates. This choice of coordinates is quite 
natural for analyzing radiation phenomena, since for hyperbolic field equations 
perturbations of the field essentially travel along characteristics. 

The most delicate part of the whole project is to gain information on the 
decay. Obviously, we cannot prepare the systems of interest and let them go to 
observe the decay of the field. In fact, at present the main interest in our whole 
subject is to obtain information on wave forms mathematically to possibly help 
identify gravitational radiation in the signals supplied by the measuring devices. 
To prepare the systems and study their evolution for an infinite time by purely 
mathematical means was out of the question at a time when not even the Cauchy 
problem local in time was completely understood. Therefore the authors above 
had to rely on their experience with a few model examples, linearized theory, 
assumptions of simplicity, physical intuition etc. to make a reasonable guess. 
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In P], Q the asymptotic behaviour was imposed by requiring that cer- 
tain metric coefficients behave in Bondi coordinates similarly as in the case of 
Minkowski space while in |3l| fall-off conditions for certain components of the 
Weyl curvature were specified. In both cases the assumptions proved consistent 
with the first few steps of a formal expansions involving the field equations. The 
authors could associate certain asymptotic data convincingly with the gravita- 
tional radiation field, and conclude from the behaviour of the Bondi- energy- 
momentum that outgoing gravitational waves carry only positive energy (|]9j, 

0). 

Soon afterwards a new idea brought into the open the basic geometric struc- 



ture underlying this analysis. In 65 , [p6| Penrose proposed to characterize 
the fall- off behaviour of gravitational fields in terms of the extensibility of the 
conformal structure through null infinity. 

To illustrate the idea we begin by describing certain conformal embeddings 
of the simply connected, conformally flat standard solutions to the vacuum field 
equations Ric[g] — Xg into the Einstein cosmos (cf. also (53]). The latter is 
given by the manifold M = M x S* 3 , endowed with the product line element 
g = ds 2 — duj 2 where duj 2 = d\ 2 +sin 2 xdv 2 denotes the standard line element 
on the 3-sphere with dcr 2 = d6 2 + sin 2 9 d <fi 2 the standard line element on 
the 2-sphere. We shall in the following be mainly interested in the conformal 
extension of Minkowski space and shall use the other two examples only to 
illustrate certain arguments. 

The case A = 0: We express the metric of Minkowski space {M ~ R 4 ,.9 = 
w^v dx^ 1 dx u ) in spherical coordinates such that g = dt 2 — (dr 2 + r 2 da 2 ) and 
define a diffeomorphism <I> of M onto the open subset <E>(M) = M' = {\s+x\ < 
I s — xl < x > 0} of M such that $~ 1 is given by 

sin s sin y 

t = 



cos s + cos x cos s + cos x 

It follows that fl 2 < I ) ~ 1 *3 = g on M' with the conformal factor = cos s + cos x 
which is positive on M' . The embedded Minkowski space acquires then a con- 
formal boundary J in M, which is defined as the subset J = d M' = J~L)J + U 

of M, where the hypersurfaces J ± = {t ± x — ±tt, < x < 
which can be thought of as being generated by endpoints of null geodesies, rep- 
resent future resp. past null infinity, the points V s1 — {x = 0, r = ±7r} represent 
future resp. past time-like infinity and the point i = {x — 7T, t = 0} repre- 
sents space-like infinity. The key point here is that the embedded Minkowski 
space together with its conformal structure and the conformal factor fl can be 
extended smoothly to the manifold with boundary M — M' U J + U J~ such 
that f2 = 0, dfl ^= on the null hypersurfaces J ± ~ K x S 2 . 

It may be of interest to note that H. Weyl already in 1923 considers the 
conformal embedding of Minkowski space into a quadric in the five dimensional 
real projective space in a relativity inspired discussion of the Mobius group and 
refers to the set which corresponds to J as to the 'sphere at infinity' (p(|). 
The fact that Weyl did not see, in spite of its beautifully geometric nature, the 
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generalization proposed in |6q] , indicates to what extent the kind of questions 
being asked and, in fact, our whole view of the subject have changed since then. 
The case A < 0: The de Sitter space is given by (M = K x S 3 ,g = dt 2 - 

cosh 2 tduj 2 ). The map M 3 (t, ti) ■% (s = arctane* - f, 0) G M, where i? 
denotes the points in S 3 , maps de Sitter space diffeomorphically onto the open 
subset M' = [— 7r/4, tt/4] x S 3 of the Einstein cosmos. It satisfies ft 2 &~ u g = g 
on M' with the positive conformal factor £1 = cos (2 s). Again the embedded 
space together with its conformal structure and the conformal factor Q can be 
extended smoothly to a manifold with boundary M = M' U J + U J~ such that 
O = 0, dd 7^ on the space-like hypersurfaces J ± = {±7r/4} x S 3 of M. 

The case A > 0. In spherical coordinates we obtain anti-de Sitter cover- 
ing space in the form (M = R. 4 ,<? = cosh 2 reft 2 — dr 2 — sinh 2 r da 2 ) . The map 

M 3 (t, r, 0, 0) (s = t, x = 2 arctan (e r ) — 5, 0, 4>) G M defines a diffeomor- 
phism of M onto the subset M' = {x < tt/2} of the Einstein cosmos which 
satisfies il 2 < E >_1 *.g = g on M' with the positive conformal factor fl = cosx- The 
embedded space together with its conformal structure and the conformal factor 
Q. can be extended smoothly to a manifold with boundary M — M' U J such 
that fl — 0, dfl 7^ on the time-like hypersurface J = {x = ~ M x 5 2 of 
M. 

The conformal extensions above are of interest already because they provide 
detailed information about the global causal relations in the three space-times 
and allow us to discuss the asymptotic behaviour of solutions to conformally 
invariant equations such as the Maxwell equations easily and with any desired 
precision. For our present discussion they are of particular interest because 
the three conformally embedded spaces and the given conformal factors provide 
solutions of the metric conformal field equations which extend, as solutions, 
through the conformal boundaries. The question whether this fact can be con- 
sidered as a model for the behaviour of more general gravitational fields leads 
us (in hindsight) to the generalization proposed in Q . 

The embedding formulae given above are distinguished by their simplicity 
but may look otherwise somewhat arbitrary. It turns out that we do not have 
much freedom. For given •& £ S 3 the curve E 9 r — > (s = 2 arctan^, i?) G M 
is a conformal geodesic on the Einstein cosmos. If x($) 7^ 0, ir, this curve inter- 
sects the conformal boundary ^ of the embedded Minkowski space for values 
T ±(^) 7^ with t_ = — r + . By specifying suitable initial data on the hypersur- 
face {t = 0}, this curve can be constructed completely in terms of Minkowski 
space and we find that the conformal factor O obtained by solving (^0|) along 
this conformal geodesic remains positive on Minkowski space and vanishes pre- 
cisely as t — ► t±, where 0(r±) 7^ 0. It follows that, apart from a coordinate 
transformation and a rescaling with a conformal factor to which is positive on 
M, the conformal extension can be constructed completely by integrating a suit- 
able congruence of conformal geodesies and an associated conformal factor on 
Minkowski space and extending the resulting formulae by continuity (cf. |4lf 
for details). 

We could have surmised the possibility for this already when we encountered 
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formula (44). It shows that for suitable initial data has precisely two zeros. 
This also raises a question. Does the fact that the formula is valid on a general 
vacuum solutions suggest that a similar procedure works under much more 
general assumptions? We are thus led again to the generalization proposed 
in g§. 

Stating in general terms what has been found above, one arrives at the 
following definition. 

Definition 3.1 A smooth space-time (M,g) is called asymptotically simple if 
there exists a smooth, oriented, time- oriented, causal space-time (M, g) and on 
M a smooth function Q, such that: 

(i) M is a manifold with boundary J , 

(ii) O > on M \J and Q = 0, dO, ^ on J, 

(Hi) there exists an embedding $ of M onto <f>(Af) = M\J such that fl 2 < I ) ~ 1 *<? 
.'/■ 

(iv) each null geodesic of (M,g) acquires two distinct endpoints on J . 

This definition, which only involves the conformal structure of (M, g) , re- 
flects the differential geometric aspects of the situations considered above. It 
implies restrictions on the global structure of the space-time. Condition (iv), for 
instance, implies that all null geodesies are complete. The set J thus represents 
a boundary at null infinity which can be thought of as being generated by ideal 
endpoints of null geodesies. Conditions (i), (ii) are somewhat redundant, since 
(ii) implies already that J is a smooth hypersurface. However, this is not the 
only purpose of (ii). Together with (Hi) it specifies how precisely $ -1 *<7 is to 
be rescaled to obtain a smooth, non-degenerate metric. The definition above 
was motivated by the following new idea (p5|). 

Penrose proposal: Far fields of isolated systems behave like asymptoti- 
cally simple space-times in the sense that they can be smoothly extended to null 
infinity, as indicated above, after suitable conformal rescalings. 

It should be noted that we have only tried to indicate the idea here. It 
should be applied with understanding and the definition above may require 
modifications in certain situations. For instance, in applying the definition to 
physically realistic situations it may be advisable to replace (iv) by some other 
completeness conditions (cf. ]50[| for an example) to allow for the possibility 
that some of the null geodesies enter black holes. The standard example here 
is the Schwarzschild-Kruskal solution, which has two asymptotic ends which 
admit smooth conformal extensions (cf. B. Schmidt, this volume). We consider 
any such changes as minor adaptations. As we shall see below, there may be 
required much more subtle ones. 

While in the earlier description of fall-off conditions coordinates play a cru- 
cial role and while even in case of asymptotically simple space-times it is often 
convenient to use Bondi coordinates, the present characterization is completely 
independent of any distinguished coordinate systems and therefore much more 
flexible than the earlier ones. But the beauty of the proposal is that it brings 
out clearly the key geometrical structure. 
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Central and critical is the requirement (iii) and the smoothness assumption. 
The smoothness properties of the rescaled field g are an expression for the fall- 
off behaviour of the physical field g. In view of the Penrose proposal definition 
( |3.1| ) acquires suddenly a much stronger interest. It is now in competition with 
another requirement and involves an extremely sharp characterization of the 
decay of fields which are governed by the quasi-linear, gauge hyperbolic field 
equations. 

The question whether the assumption that the field equations hold near J 
is consistent with the smoothness and conditions (i) - {iii), or whether the 
smoothness assumptions made in definition need to be relaxed has been 
and, as seen below, still is a matter of debate. The question is not whether C°° 
should be replaced by C k for some large k. The question is whether solutions to 
the field equations admit conformal extensions of class C k where k can be chosen 
large enough to make the proposal meaningful and, related to this, whether the 
conformal Weyl tensor tends to zero on the conformal boundary. 

Penrose addresses this problem in J66| . For later comparison we recall the 
assumptions and give (a version of) the main argument, using the spin frame 
formalism. Consider a smooth solution to the metric conformal field equations 
with A = on the physical manifold M and assume that it admits a conformal 
extension to a manifold M = M U J + of class C 4 such that g and CI extend to 
fields in C 3 (M). The conformal Weyl spinor ^> a b c d is then in C 1 (Af). 

By definition, CI = 0, dCl ^ and, as observed earlier, V aa <fiV aa CI = 
on J + . We assume that there is a space-like surface E in J + which intersects 
each null generator of J + exactly once. In a neighbourhood of E there are 
then C 3 spinor fields o a , i a with i a o a = 1 such that V QQ ' Cl = -o a d a ' on J+ 
and the complex vector field m aa = o a I a is tangent to E. Note that the 
set of such spin frames defines a reduction [7(E) of the bundle of normalized 
spin bases over E to the subgroup [7(1) of SL(2,C), where e*3 g [7(1) acts 
on the spin frames by (o,i) — * (e 1 ? o,e _I 5t). A 2 : 1 bundle morphism maps 
[7(E) onto the bundle + (E) of oriented frames which are orthonormal for 
the induced metric h on E such that we can write m = l/y/2 (n + ir-z) with 
r = {va}a=i,2 S + (E). Denote by u the generator of the structure group 
SO(2) of + (E) with components u A B = S A 1 S' 2 B - S A 2 S 1 B and by Z the 
vertical vector field on + (E) generated by u. 

Equation ( |3^ ) reduces on J + to V oa ' V&6' CI — se a bE a >b>- Since o a V aa ' is 
an interior operator on J + it follows that o a V aa '(o( ) Of ) /) = —sohta'V on J + 
from which we obtain by suitable contractions i a V m o a + r° V TO d a / whence 
7 = g(S7 m m,fh) = — 2i a V m o a . Given 7, the connection T> of h is known, 
because T> m m = —'ym, T> m ffi — jfh. 



With (f>abcd = CI ^abed the Bianchi equation (36) takes on M the form 



v Q a , n v abcd - n v a „/ y abcd = 0. (47) 

(We note that in the case where ^ A = — 3V//'f2V77 on J this implies 
at each point of J immediately that = — 3 V ea <f2 V aa a bcd = A ^ e bcd-) It 
follows that $ abed — cto a Ob o c Od on J + with a — ^> a bcd t a i b i c i d ■ If in addition 
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it is assumed that 



VLV ee ,V a a'^abcd^O at J+, (48) 

then equation ( p7| ) implies, after taking a covariant derivative, symmetrizing in 
a' and ef, and performing some spinor algebra, the equation o e V ee ' ^ a bcd = 
on J + . With the special form of the conformal Weyl spinor on J + it follows 
after contraction with tr 

V m a-2ja = on J+. (49) 

By its definition the C 1 function a can be considered as a 0-form on + (E) of 
spin weight —2, i.e. 

Za = -2ia on 0+(E). (50) 

Equation ( |49| ) says that the covariant differential of a on + (E) (cf. plj) 
vanishes on vector fields which project onto m. It can thus be expressed in the 
form 

(Hi + iH 2 )a = on 0+(E), (51) 

where for A = 1, 2 we denote by Ha the horizontal vector field on + (E) which 
projects at the point r £ + (E) onto the tangent vector ta of E. 

Suppose now that E is diffeomorphic to the two-sphere S 2 . Because S 2 carries 
only one conformal structure the conformal scaling of g can be chosen near E 
such that (E, h) can be identified isometrically with (S 2 , — da 2 ). This allows 
us to identify + (E) with SO(3), the isometry group of S 2 , such that there 
are left invariant vector fields Y t , i = 1,2,3, on SO(3) which satisfy under this 
identification Y\ = H2, Y% — H±, Y3 = Z and [li,Yj] = E^yj, where e^fc 



is totally anti-symmetric with 6123 = 1. Equations (p0|), (51) then imply the 
distributional equation 

-\{Y? + YZ+Y 2 )a = a on 0+(E). 

Since the operator on the left hand side is the Casimir operator with respect to 
the Killing form on 50(3), the solution space of this equation is well known (cf. 
[[76)). It is spanned by the matrix elements of an irreducible three-dimensional 
unitary representation of the group ST/ (2) (in the notation of IsS] the functions 
T2 1 j), which descends to SO(3). Because the only function in this space which 
satisfies ( |50"1 ) is a = 0, it follows that the Weyl spinor field vanishes on E. Thus, 
4> a bcd extends in a continuous fashion to and, as a consequence, satisfies 
Sachs' peeling property (^36j, [[73], @)- 

Remarkable as it is that such a conclusion can be drawn from the special 
conformal covariance of the Bianchi equations, the fact that we are dealing 
with a spin-2 field, and the topological assumption on E, there remains the 
question whether it follows from the nature of the long time evolution by the 
field equations that the smoothness assumptions made above and (Q) can be 
satisfied and whether this is in fact the case for a sufficiently large class of 
solutions. 
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There are two different issues here. We start with a discussion of the situation 
where the vacuum field equations with a cosmological constant A are satisfied 
locally near a point p of J and where fl and g are smooth. As seen earlier, the 
conformal field equations then imply V p fiV p fi = —A/3. In particular, J is a 
null hypersurface if A = 0. In that case we can expect J to have two components 
J ± and p £ J + . To get some idea of the class of solutions which are compatible 
with conditions (i) - {in), we consider the asymptotic characteristic initial value 
problem for the metric conformal field equations where data are prescribed on 
two intersecting hypersurfaces which are to assume, by the arrangement of the 
initial data on them, the meaning of an outgoing null hypersurface Af which 
meets J + in a space-like 2-surface S and the part of J + in the past of S. 
This is essentially the problem underlying the formal expansions studied in || , 

till- & 

All solutions for which the rescaled conformal Weyl tensor d l ju has a smooth 
limit on can be characterized in the analytical as well as in the C°° case 
(p^j, ]57t) in the sense that for suitably given data the existence of a unique 
local solution near £ can be shown. The freedom to prescribe data on and 

, which are in each case two components of d 1 jki referred to as the null 
data, turns out to be the same as in the standard characteristic initial value 
problem for Einstein's vacuum field equations (corresponding to the special ge- 
ometric nature of the hypersurface , there is a slight difference with the 
latter problem concerning the data which may be prescribed in addition to the 
null data on the intersection of the two null hypersurfaces). If, as above, the 
tensor d l jki is represented by the symmetric spinor (f> a bcd, then (j) a bcdt- a t h f t d 
and 4>abcd o a o b o c o d denote the null data on Af and J + respectively. It is quite 
a remarkable property of the conformal field equations that they define the con- 
formal boundary {VL = 0} in terms of one of their unknowns without restricting 
the freedom to specify smooth data there. 

The null data (f> a bcd o a o b o c o d on J + have a natural interpretation as the 
outgoing radiation field. Further important physical concepts can be associated 
with the hypersurface or subsets of it and the questions of interpretation 
have been extensively analyzed (cf. also (5), |4j|). All indications suggest that 
we are on the right track. 

The situation is similar in the de Sitter- type case A < 0. Here null infinity 
is space-like. We ask whether we can characterize the solutions which admit 
a smooth conformal extension at past null infinity J~ under the (natural) as- 
sumption that this hypersurface is orientable and compact as in the case of de 
Sitter space and the further assumption that d l jki extends smoothly to J~ . 
Again, for suitably given data on J~ there can be shown the existence of a 
unique solution to the conformal field equations. The freedom to specify data 
on J~ is essentially the same as in a standard Cauchy problem with the ex- 
ception that the mean intrinsic curvature of the initial hypersurface cannot be 
prescribed (]29[|). This appears natural, because J~ is a very particular hyper- 
surface with respect to the solution. Another peculiar feature consists in the 
fact that the Hamiltonian constraint becomes trivial on J~ . 

The situation where A > 0, in which case null infinity is time-like, differs 
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in various aspects from the other cases. In this case the assumption of the 
smoothness of this limit has a different status than in the other cases. It appears 
that in general there would not even exist complete null geodesies if the field 
d l jki tended to blow up sufficiently fast 'far out'. However, there is available 
a detailed discussion of the existence of solutions for which d l jki has a smooth 
limit to J (Q). 

The solutions studied in these three situations are semi-global in the sense 
that in terms of the physical metric their domains are of infinite extent and 
even comprise asymptotic regions. However, this is not sufficient to answer the 
question posed above. Here comes in the second issue concerning the compati- 
bility of asymptotic simplicity with the field equations. What we need to control 
is the behaviour of the solutions as they evolve towards null infinity. Thus, we 
cannot assume anything about the field d % jki (except in the anti-de Sitter-type 
solutions where boundary data must be given at space-like and null infinity 
to characterize the solutions uniquely) and it is not clear a priori whether it 
remains bounded near null infinity. We need to find out whether it behaves, 
together with the other unknowns in the conformal field equations, such that a 
conformal extension (M, g) of sufficient smoothness can be constructed. 

There exists a model case which illustrates the type of result we would like to 
have. In the case of the de Sitter space considered above (A = —3) we consider 
the hypersurface S = {t = 0}. The Cauchy data induced by de Sitter space on 
5* are given by the standard metric on S 3 and a vanishing extrinsic curvature. 
It turns out that data for the equation Kic[g] = Xg which are sufficiently close 
to the de Sitter data evolve into solutions which are asymptotically simple in the 
sense of definition ( |3 . 1[ ) , or, in other words, asymptotic simplicity is non-linearly 
stable p2[). Here d i ju attains a smooth limit at the conformal boundary as 
a consequence of the evolution process and the assumptions on the initial data. 
The fact that the Schwarzschild-de Sitter solution has time slices diffeomorphic 
to S 2 x S 1 but develops horizons, singularities, and only patches of a smooth 
boundary shows that neither the compactness of the initial hypersurface nor the 
sign of the cosmological constant is sufficient for this behaviour. Remarkably, it 
has been shown recently that a solution with a compact Cauchy hypersurface can 
admit smooth conformal boundaries in the past as well as in the future only if 
the Cauchy hypersurface satisfies the topological restriction that its fundamental 
group be finite (Q], cf. G. Galloway, this volume, for further results concerning 
de Sitter- type solutions). 

An analogous stability result in the Minkowski-type case A = is still lacking. 
We shall see that a result of this simplicity cannot be expected. 

4 Asymptotic behaviour of vacuum fields with 
vanishing cosmological constant 

For physical applications one would like to know how the sources interact with 
each other and with the gravitational field and to what extent properties of 
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the far field reflect this interaction. For a first mathematical analysis of the 
general behaviour of the far fields details of the sources are irrelevant and it is 
natural to ignore them in a first study. Thus we are led to ask whether there 
exist other solutions to the Einstein equations Ric[<7] = than Minkowski space 



which satisfy the requirements of definition (3.1 



Two different types of initial value problems have been studied in this con- 
text. The most important one is the standard Cauchy problem where Cauchy 
data are prescribed on a space-like hypersurface S with one or several asymp- 
totically Euclidean ends. If there are several such ends we can expect black 
holes to develop and condition (iv) to be violated. Therefore, and because it 
is the evolution near an end we want to focus on, we will concentrate on the 
case S ~ K 3 with one end. Standard examples are the hypersurface {t = 0} in 
Minkowski space with the induced data or perturbations thereof. 

The other type of problem is the hyperboloidal initial value problem where 
data are prescribed on a hypersurface H with boundary E, which is thought of as 
a hypersurface in an asymptotically simple space-time which extends to future 
null infinity 3 + , intersects J + in E, and is space-like everywhere. The boundary 
E can have several components, examples exist in the Schwarzschild-Kruskal 
space-time (cf. B. Schmidt, this volume). Again we will concentrate on the case 
where only one asymptotic region exists such that i/~{x€lR 3 ||a;|<l} and 
E ~ S 2 . The standard example of such a hypersurface is given by the extension 
of the space-like hyperboloid {t 2 — \x\ 2 = l,t > 0} in Minkowski space to J + . Its 
interior inherits from Minkowski space a metric of constant negative curvature 
which can be rescaled to become Euclidean. 

While the hyperboloidal initial value problem has been introduced initially 
as an intermediate step towards a treatment of the standard Cauchy problem 
( pff] ) , it has gained in recent years an interest in its own as a basis for numerical 
calculations (cf. J. Frauendiener and S. Husa, this volume). It is intrinsically 
non-time-symmetric. The future evolution will extend under suitable assump- 
tions to null infinity while the past time evolution develops a Cauchy horizon 
which corresponds to a part TV of an outgoing wave front. The domain covered 
by its evolution coincides with that which can be determined by the characteris- 
tic initial value problem where data are prescribed on N (if this set is sufficiently 
smooth) . 



4.1 The hyperboloidal initial value problem 

Since, for reasons which will become clear later, the hyperboloidal initial value 
problem is much easier to analyze in the context of asymptotic simplicity than 
the standard Cauchy problem, it has been analyzed much earlier and much 
more is known about it. If one wants to construct its solutions by solving the 
conformal field equations, the first task is to provide conformal initial data. In 
the analyses available so far, in which a non-vanishing constant mean extrinsic 
curvature is assumed, certain data can be prescribed on H arbitrarily, certain 
data are obtained by solving elliptic equations which degenerate on E, and the 
remaining data are obtained by differentiation and algebraic operations which 
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involve divisions by the conformal factor, which vanishes on E. We give an 
outline of the results needed in the following discussion (for further information 
cf. L. Andersson, this volume). 

The subtleties arise at X. In CJ, [|| free data are prescribed which extend 
smoothly to S. It turns out that (a) in general the remaining data have expan- 
sions in terms of powers of r and logr, where r denotes the conformal distance 
from E (i.e. they admit polylogarithmic expansions at E), while (b) all data will 
be smooth if the free data satisfy certain regularity conditions which require that 
a few functions derived from the free data vanish on E. A much more general 
class of data is considered in In particular, (c) free data are prescribed which 
admit polylogarithmic expansions at E and it is shown that the remaining data 
also admit such expansions. 

The conformal field equations have been used to study the evolution of some 
of these data. In the case (b) the data evolve into solutions which satisfy the 
first three conditions of definition (|3.l| ) and which admit in fact a conformally 
regular point i + at time-like infinity, analogous to the point i + in the case of 
Minkowski space, if the data are sufficiently close to Minkowskian hyperboloidal 
data (|2(|), |50|, |3^|). Recently also the existence of solutions to the conformal 
field equations for data of the type (a) has been shown on a manifold which 
contains a boundary hypersurface J' on which the unknown of the conformal 
field equations vanishes. The solution is 'rough' at J' and it is expected that it 
admits a polyhomogeneous expansion there (pc|], cf. also iQ). The evolution 
of data of type (c) has not been studied yet. 

While the basic idea underlying the Penrose proposal, namely to character- 
ize the asymptotic behaviour of solutions in terms of conformal rescalings, is 
still being used in this study of asymptotically rough data by employing the 
conformal field equations, this situation certainly forbids any straightforward 
application of definition (3.1). The result seems to indicate that the smooth- 
ness assumption in definition ([O]) may be too strong. However, before we can 
arrive at a conclusion here we need to answer the following question: of which 
physical relevance are data of type (a) and (c), are they needed to model the 
systems of interest and which are the 'systems of interest'? 

A clear answer could be given, if it could be shown that standard Cauchy 
data by necessity develop into solutions containing hyperboloidal data of type 
(a) or (c). The results for data of type (b) tell us that the decision about the 
smoothness of the conformal structure near null infinity is essentially made near 
space-like infinity. Thus we need to analyze the behaviour of the solution to the 
standard Cauchy problem there. 



4.2 On the existence of asymptotically simple vacuum so- 
lutions 

We shall see below that the general analysis of the conformal structure of solu- 
tions to the Einstein equations near space-like infinity is quite complicated. On 
the other hand, static or, more generally, stationary asymptotically flat vacuum 
solutions are known to admit smooth conformal extensions through null infinity 
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(cf. p8[). Thus, to avoid the difficulties discussed below one could try to con- 
struct asymptotically flat standard Cauchy data (on R 3 say) which agree with 
stationary data near space-like infinity. The time evolution of these data would 
then be stationary near space-like infinity and clearly contain smooth hyper- 
boloidal hypersurfaces. One could then hope to combine this with the results 
of |30j to show the existence of solutions to the Einstein equations which satisfy 
all conditions of definition ( |3.1[ ). 

The first attempt at this was made by Cutler and Wald who succeeded 
in constructing smooth initial data for the Einstein-Maxwell equations with 
arbitrarily small ADM mass which agree with Schwarzschild data near space- 
like infinity (0). Using the results of J32|, they were able to show for the 
first time the existence of smooth solutions to the Einstein-Maxwell equations 
which satisfy the conditions of asymptotic simplicity. The vacuum case remained 
difficult. 

More recently Corvino developed a gluing construction by which an arbitrary 
smooth, time-symmetric, asymptotically flat solution (R 3 , h) to the constraints 
can be deformed on a certain annulus {x S R 3 |0 < R < \x\ < 2R}, such 
that it extends with a given finite differentiability to Schwarzschild data out- 
side 2 R while remaining unchanged inside the radius R and still satisfying the 
constraints everywhere (|l6||). It turned out difficult, however, to keep the ra- 
dius R fixed or to control its growth sufficiently if the construction is applied to 
families of initial data for which the ADM mass m is going to zero. This makes 
it impossible to obtain the desired result by an immediate application of p0[ |. 

Quite recently Chrusciel and Delay obtained a similar result based on a 
modification of Corvino's technique in which the radius R as well as the differ- 
entiability class C k may be kept fixed while m — > if they base their construction 
on initial data (R 3 , h) such that in the standard coordinates on R 3 the parity 
condition h(x) = h(—x) is satisfied by the 3-metric (p5[|). In this case it is 
easy to construct in a systematic way hyperboloidal hypersurfaces in the time 
evolution of the data which contain the part {x £ R 3 | |a;| < 2 R} of the initial 
Cauchy hypersurface as a subset and which extend into to future of the initial 
hypersurface in the Schwarzschild part of the evolution. It can be arranged that 
the induced data approach Minkowskian hyperboloidal data if h approaches the 
Euclidean metric on R 3 (cf. [|d| for a construction based on conformal Gauss 



systems and 38 for the discussion of a smooth limit of Schwarzschild solu- 
tions to the Minkowski solution near space- like infinity). Since this can be done 
with sufficient smoothness, the results of |3(J and the fact that the data are 
time-symmetric then imply the existence of non-trivial solutions to the vacuum 
field equations which satisfy the conditions of definition (3.1) with a finite dif- 
ferentiability of the rescaled metric which can be chosen to be large. For these 
solutions the rescaled conformal Weyl tensor is differentiate on the conformally 
extended space and there exist conformal extensions containing regular points i 
which represent future and past time-like infinity. This resolves a question which 
remained open and controversial for forty years. 

While it is reassuring to see that we are not discussing a class of solutions 
which consists only of Minkowski space, it is clear that the result does not 
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answer all our questions. The solutions obtained here are very special. Being 
Schwarzschildean near space-like infinity, their NP constants vanish ( jl2| , ]62|) . 
As a consequence the rescaled conformal Weyl tensor d l jkl vanishes at (]46[, 
[|63|). However, from the numerical calculation of solutions with a regular i + 
we know that in general d l ju does not vanish at i + (|54 ). This suggests that 
modifying initial data sets to become Schwarzschild near infinity considerably 
reduces their radiation content. 

Hopefully, there will be developed generalizations of the method which will 
allow us to glue other stationary data to given data such that we prescribe 
multipolcs moments and non- vanishing NP constants. Then the argument above 
will not be available any longer and the data will gain in interest in the modeling 
of physical systems. 

But this will still not help us answer the decisive question: what degree of 
restriction is implied by 'asymptotic simplicity' in the class of all ('reasonably' 
defined) asymptotically flat solutions. 



4.3 The regular finite Cauchy problem 

Under the conformal embedding of Minkowski space into the Einstein cosmos 
described above, the Cauchy hypersurface 5" = {t = 0} ~ R 3 with its interior 
Euclidean metric e is mapped conformally into the unit 3-sphere, S — ► S = {s = 
0} = SU {i}, with its standard metric duj 2 and a point i with coordinates s = 0, 
X = 7r which represents space- like infinity for (S 1 , e). If more general data h, x 
satisfying the vacuum constraints are prescribed on S which are asymptotically 
Euclidean in the sense that h suitably approaches the metric e near space- 
like infinity it appears reasonable to conformally compactify («S, h) by suitably 
rescaling h — > h = il 2 h and adding a point to S as above. Depending on h. 
there may not exist coordinates x a on S near i and a conformal factor f2 such 
that the metric coefficients h a /3 will be smooth at i. But even if the data h, x 
are chosen such that h extends smoothly to i, in fact, even if x — 0, it turns out 
that in coordinates adapted to h the rescaled conformal Weyl tensor behaves as 

d^x P = 0(r- 3 ) as r -> 0, (52) 

where r denotes the /i-distance from i, unless the ADM mass of the data van- 
ishes, i.e. h is flat. 

To let things not look utterly hopeless, we note that d M v \ p is not an arbi- 
trary tensor with Weyl symmetry satisfying (p2[), that its structure near i is 
restricted by the assumption that h and x be asymptotically Euclidean, satisfy 
the constraints, and that v \ p is given by a specific formula. We refer to |38j 
for explicit general expressions. 

In the case of Minkowski space the hypersurfaces J ± , which are swept out 
by the future resp. past directed null geodesies emanating from the point i° 
(naturally identified with i G S), are characteristics for the conformal field 
equations in the conformally extended space-time. If we assume for a moment 
that a similar picture can be established for the space-time which develops from 
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non-trivial data h, x, then ( j5^ ) makes us wonder why this strong singularity at 
i does not spread along the characteristics J ± and destroy the smoothness of 
the conformal structure there in the first place ! If such a contradiction does 
not arise, it can only be due to (a) specific features of the field equations and/or 
(b) specific requirements on the data which are analogous to the regularity 
conditions observed in the hyperboloidal initial value problem. 

If features as referred to in (a) exist, they must lie beyond conformal reg- 
ularity. If requirements such as those referred to in (b) exist we may wonder 
what kind of conditions could, in a situation where the conformal extension is 
only given by a point, replace the vanishing of certain functions on a 2-surface 
as required in the hyperboloidal problem. Moreover, if one tries to analyze an 
initial value problem for the conformal field equations for data as singular as 
( [52| ) the choice of a gauge becomes a delicate problem and there is little lee-way 
for rough estimates. 

Somewhat unexpectedly it turns out that for suitably chosen asymptotically 
flat initial data the standard Cauchy problem can be reformulated as a regular 
finite initial value problem for the general conformal field equations with data 
given on a compact manifold with boundary. Moreover, the underlying setting 
also allows us to analyze to some extent the questions raised above. 

The underdetermined quasi-linear system of constraints imposes very weak 
conditions on the data. In particular, it requires very little with respect to the 
smoothness of the data at space-like infinity. Thus it is clear that we have to 
make assumptions concerning the precise behaviour of the data near i. In J38|, 
were the setting has been developed and worked out in detail as described below, 
it was assumed that the data are time-symmetric, that h extends smoothly to i, 
and that h is in fact analytic in normal coordinates at i. It follows in particular 
that x = an d that the conformal factor f2 which relates h to the physical 
metric h = fT 2 h satisfies Q € C 2 {S) UC°°(S\ {i}) with Q = 0, cM = 0, 
Hess f2 ~ h at i, but is not smooth at i unless the data have vanishing ADM 
mass. As a consequence we still have (|52|). It should be emphasized that our 
assumptions have been made for convenience because the structure of the data 
near space- like infinity needs to be known in all details (cf. the remarks below) . 
In particular, the assumption of analyticity can be dropped and replaced by a 
C°° condition. 

4.3.1 The cylinder at space-like infinity 

In the following we shall describe the main considerations leading to the for- 
mulation of the regular finite initial value problem near space-like infinity. The 
general conformal field equations will be used here for several reasons. We want 
to obtain a formulation which is leading to statements which depend as far as 
possible on the conformal structure itself and as little as possible on conditions 
'put in by hand'. The nature of the singularity at i forces us to carefully dis- 
tinguish between the dependence of the fields on radial and angular directions 
at i. The frame field in the formalism is readily adapted to this situation. 
Much more importantly, we want to avoid any gauge which depends on implicit 
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wave equations, because these might transport non-intrinsic singularities along 
null infinity. An example for this occurs if the conformal gauge in the metric 
conformal field equations is removed by prescribing the Ricci scalar. Clearly, 
space-like infinity is not met by time-like conformal geodesies starting at points 
of S = S \ {i} and conformal Gauss systems should thus be unaffected by the 
singularity at i. Finally, the location of the set where the conformal factor 
vanishes will be known a priori. 

In |3^| the general conformal field equations have been used in the spin frame 
formalism, because this leads to various algebraic simplifications. Here we will 
discuss the standard frame formalism, suppressing many details for which we 
refer the reader to Q . 

Assume a fixed choice of the conformal scaling for the initial data. Choose 
a fixed oriented ft,-orthonormal frame e a , a = 1, 2, 3 at i. Any other such frames 
at i is then obtained in the form e a (s) = s c a e c with s = (s c a ) € 50(3). For 
given s, we distinguish es(s) as the radial vector at i and assume e a {s) to be 
parallelly transported with the Levi-Civita connection of h along the geodesies 
with tangent vector es(s) at i. We denote by p the affine parameter on these 
geodesies which vanishes at i and denote the frame obtained from e a (s) at 
the value p by e a {p 1 s). It will be assumed that \p\ < a, where a is chosen 
such that the metric ball B centered at i with radius a is a convex normal 
neighbourhood for h. Then ] — a, a[xS0(3) 3 (p, s) — ► e a (p, s) g SO(S) defines 
a smooth embedding of a 4-dimensional manifold into the bundle 50(5) of 
oriented orthonormal frames over 5. 

We denote by B the image of the set [0, a[x 50(3), by 1° its boundary 
{p = 0} ~ 50(3), and by tt the restriction to B of the projection of 50(5) onto 
5. For given s G 50(3) the vectors e3(ss') are identical for s' G 50(2), the 
subgroup of s = (s c a ) G 50(3) for which s c 3 e c = e^, and all these vectors are 
parallelly transported along the same geodesic. It follows that the set B inherits 
from 50(5) an action of the subgroup 50(2) which implies a factorization 

B^B 1 = 5/50(2) ^ B of the map tt such that tt" maps the set ~ 5 2 

onto i while it implies a diffeomorphism of B' \ ) onto the punctured ball 
B \ {i} which we could use to identify these sets. 

It turns out more convenient, however, to pull the initial data on B via tt 
back to .B, use this set as the initial manifold and p and s as 'coordinates' on 
it. Questions of smoothness are then easily discussed and certain expressions 
simplify considerably. Since all fields have a well defined transformation be- 
haviour (spin weight) under the action of 50(2) and this action commutes with 
the evolution defined by the propagation equations, no problem will arise from 
this. In ]3£| ] it is shown how to prescribe on the part of B where p>Oa frame 
field X, c a (p,s), a — 1,2,3, where X is generated by the action of 50(2), the 
vector fields c a (p,s) satisfy T(ir)c a (p, s) = e a (p, s), and a frame formalism can 
be applied. In the coordinates p and s on B, which are essentially spherical 
coordinates, some frame vector fields and certain data derived from them are 
singular at p = 0. 

With a suitable choice of the initial data for the conformal geodesies the 
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conformal factor and the 1-form df. take the form 

6 = -fl-r 2 4V d = -2r^, d a = -(e Q (0)), (53) 
where = 2 fi |/i Q/3 ft Q fi ,/5.L! /2 and 

k is a function on B which we are free to 
choose. Any function /* in ( p3|) is considered to be constant along the conformal 
geodesies and to agree with its given initial value /. In the conformal extension 
of Minkowski space in which space-like infinity is represented by the point i° we 
would have k = O(l) as p — > 0. For our purpose it will be much more useful to 
choose K in the form k = pn' with some smooth scalar (i.e. S*0(2)-invariant) 
function n' on B with k' = 1 on 1° such that 0| T= o = — = O(p) as p — > 0. 

In this particular conformal gauge all data extend smoothly to 1° C B and the 
coefficient functions 0, dk in the general conformal field equations are smooth 
as well. 

We are free to choose k! . For definiteness we set k = w such that = 
for t — ±1. If the resulting coordinates remain regular long enough, future 
and past null infinity will be represented near space-like infinity by the sets 
J± ={t = ±1}, while space- like infinity will now be represented by the cylinder 
at space-like infinity which is given by the set I = {p = 0, |r| < 1}. It 'touches' 
null infinity at the sets I = {p = 0, r = ±1}. What in the beginning was a 
Cauchy problem now looks like an initial boundary value problem. It turns out, 
however, that the boundary / is of a very special nature. 

If we write w = (e M F; J : Ljk) and the unknown ( p8[ ) in the form u = (w, z) 
with z — (d l jki), the symmetric hyperbolic reduced equations take the form 

A> i (w)d ll z = H(w)z, d T w = F{w,z), (54) 

(where we suppress the dependence of F on the conformal factor and the 1- 
form dk)- Extending the data on B as well as the evolution equations (preserving 
symmetric hyperbolicity) smoothly into the domain p < 0, standard results 
imply the existence of a smooth local solution in a neighbourhood of the initial 
hypersurface which contains a piece of /. Fhe fact that = on / then allows 
us to conclude from the details of the equations above that 

A p = on I. (55) 

This implies that / is a total characteristic in the sense that the complete system 
( |54"| ) reduces on / to a symmetric hyperbolic system of inner equations on /. It 
follows that u is determined uniquely by its values on 7°, that the solution 
obtained for p > does not depend on the chosen extension of the data through 
7°, and that u extends smoothly to I (which can not be taken for granted in 
general initial boundary value problems, cf. Q, Q). 

The construction of I can be considered as attaching a boundary at space- 
like infinity to the physical space-time. However, we wish to emphasize that 
not the definition of a boundary but the analysis of the field equations is our 
main interest. We do not impose any conditions on the time evolution of the 
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solutions except the field equations. The construction is defined in terms of 
general properties of conformal geometry, the conformal field equations, and 
the choice of the initial data for the gauge. 

While everything extends smoothly to I near 1°, it should be observed that 
/ is not a geometrical entity but that it is attained as a certain limit of geometric 
objects, namely the conformal geodesies defining our gauge. It is not a time- 
like hypersurface and some of the metric coefficients do in fact diverge at I 
in the given coordinates. This is consistent with the smoothness of the frame 
coefficients and the vanishing of the matrix A p on /. Since the wave equation for 
the rescaled conformal Weyl tensor would degenerate at /, it turns out important 
here that we extract the hyperbolic system from the first order Bianchi equation. 
It is remarkable that the blow-up procedure i° — > /, the coordinate p which is 
adapted to the conformal scaling defined on S by f2, and the conformal gauge 
defined by O conspire to produce a regular finite problem. By a simple rcscaling 
(different choice of k) and a corresponding redefinition of the coordinates the 
representation of space-like infinity in terms of the cylinder I can be converted 
into the conventional representation of space-like infinity by a point i° (]3§||). 

4.3.2 The transport equations on the cylinder at space-like infinity 

The procedure described above, by which the point i° is replaced by the cylinder 
/, does not only give us regular data and equations but also leads to an unfolding 
of the evolution process near space-like infinity which allows us to analyze the 
process there at arbitrary order and in all details. Taking formal derivatives 
of the evolution equations with respect to p, restricting to the cylinder, and 
observing (p35|), one obtains transport equations on / for all the functions u p — 
d p u\i, p — 0,1,..., on the cylinder. Expanding u p = u p (t,s) in terms of a 
certain functions system Ti m 1 k (where m = 0,1,2,... and i,k = 0, . . . , 2 m), 
on SO(3) (resp. on SU{2) in the case of the spin frame formalism used in ]3^|), 
these equations reduce to a hierarchy of ODEs along the curves ] — 1, 1 [5 t — > 
(r, p = 0, s) 6 /, where the differential operator depends on p but not on the 
data while the right hand side of the equation for u p+1 depends quadratically 
on w°, . . . , u p and on the coordinates. 

The calculation of u on I allows us to determine A p (u) on /, which gives in 
particular 

A T = diag(l + t, 1, 1,1,1 — r) on I. 

While the second of equations (Q) is an ODE which is regular at I ± , the matrix 
A T occurring in the first of equations ( f54| ) degenerates precisely on I ± . All the 
remaining open questions concerning the regular finite initial value problem near 
space-like infinity are related to this fact. 

This degeneracy (which is removed at the next order in p) has important 
consequences. It turns out that in general the quantities u p , in particular the ex- 
pansion coefficients of the rescaled conformal Weyl tensor, contain terms which 
behave as (1 — r) k log-'(l — r) at I + . In the expressions which so far have been 
calculated explicitly only terms with the powers j = 0A occur, but due to the 
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non-linearity of the equations other powers should arise as well. Because of 
the hyperbolicity of the reduced equations these singularities can be expected 
to spread along 3 . The powers k and j observed in the explicit calculations 
then suggest that in general not all components of <f ji m will be bounded at 
null infinity. In fact, in the spinor notation used in (f47|), ( fl8| ) it appears that 
in general the rescaled Weyl spinor field will have an unbounded component 
<l>abcd t a ^ t c t d which behaves as log(l — r) at J + , where r — > 1. 

This suggests that in general the conformal Weyl tensor will vanish on J + 
but will neither be differentiable there nor satisfy the assumption (pl8|). 

Some care is required here, because the hyperbolicity does not extend to J*. 
However, if the whole setting developed in |Q is linearized at Minkowski space, 
a complete analysis of the situation can be given which confirms the conclusion 
(cf. J. Valicnte-Kroon, this volume). 

Due to the geometric nature and the conformal invariance of our gauge 
conditions we are discussing here intrinsic features of the solutions. If the null 
infinities of the solutions would admit smooth extensions of the conformal struc- 
tures, then the conformal geodesies, the conformal factor, and the frame would 
extend smoothly through null infinity. Also the components of the rescaled 
conformal Weyl tensor in that frame would then extend as smooth functions of 
r through null infinity. Thus the observations above suggest that even under 
assumptions on the conformal data which cannot be improved from the point of 
view of smoothness there will in general occur logarithmic singularities at null 
infinity due to the evolution process. 

4.3.3 A conjecture 

The class of logarithmic singularities determined in |38| have been found by 
solving the transport equations on / explicitly for certain components u' p of the 
functions u p , p = 0, 1, 2, . . ., in the decomposition with respect to the functions 
Tim 1 k- Inspecting the initial data u ,P _ for these quantities on 1° and relating 
them to the free datum h on S, we find that the logarithmic terms we find for 
general data in our class do not occur in u' p for p < + 2 if and only if h 
satisfies the regularity condition 

S(D il ,...,D iq B jk (i)) = for q = 0, 1, . . . , g*, (56) 

where D denotes the connection and Bjk — \ tj %l Di Lik the Cotton tensor of 
h while S denotes the operation of taking the symmetric trace-free part of a 
tensor. 

Consistent with the fact that the free initial datum is given by the confor- 
mal structure of h, condition ( |56| ) is conformally invariant for given q* > 0. 
The condition has been found earlier in the quite different, though related, con- 
text of vanishing ADM mass (plj). Although the derivation is much easier in 
that case and related to a geometric picture, we are still lacking a geometrical 
interpretation of condition (f56|). 

Since the condition is conformally invariant and it is easy to construct data 
which are conformally flat near space-like infinity, there exist many non-trivial 



38 



data satisfying the condition. With the techniques of |15 and such data 
can even be constructed so that they agree with prescribed data on a compact 
set. Moreover, there exists a large class of data which satisfy the condition in a 
non-trivial way (p8[). In fact, all data which are conformal to static data near 
space- like infinity satisfy (|56|) for all q* (f3jj). Since static data are analytic 



near space-like infinity (|6|, p9|), it follows that condition (|56]) can be satisfied 
with Cotton tensors which do not vanish near space-like infinity (p8fl). This is 
consistent with the fact that the static solutions represent the largest class of 
time-symmetric solutions which are known to admit near space-like infinity a 
smooth structure at null infinity. It would be interesting to know whether the 
results of [|l5) and Jl6| can be generalized to glue static ends to given data. 

For various reasons we expect condition ( |56|) to be of a much wider signifi- 
cance. 

Conjecture: There exists an integer fc* > such that for given k > k* 
the time evolution of an asymptotically Euclidean, time- symmetric, conformally 
smooth initial data set admits a conformal extension to null infinity of class C k 
near space-like infinity, if the Cotton tensor satisfies condition (^6|) for a certain 
integer = q*(k). 

If this conjecture can be shown to be correct, solutions for data which sat- 
isfy ( p6| ) for all q* will contain smooth hyperboloidal hypersurfaces. We expect 
that the corresponding hyperboloidal initial data approach Minkowskian hy- 
perboloidal data in a continuous way if the Cauchy data approach Minkowski 
data. This then can be used in conjunction with the results of pp[ | about the 
hyperboloidal Cauchy problem to show the existence of solutions satisfying all 



the requirements of definition (3.1). Because of the degeneracy of the equations 
at I ± the proof of the conjecture is not a standard problem and requires new 
ideas. However, we can split off an important and interesting subproblem. 

Sub conjecture: // condition ( |5q ) holds for a given q* > 0, the functions 
uP , P < Q* + 2, extend smoothly to 7*. 

In the formalism employed in | f38| the validity of the statement above can in 
principle be checked by a recursive calculation of the functions vP . In fact, such 
a calculation shows that the u p are smooth at I for p < 3 if the corresponding 
regularity conditions hold (f38), |4^| ). The general proof essentially reduces to 
controlling the algebraic structure of the terms on the right hand side of the 
ODEs governing the transport of the u p , which, with increasing order, becomes 
quite complicated (cf. the discussion in j42| , J43|). It is most desirable to deter- 
mine terms of higher order by an explicit calculation, possibly with an algebraic 
computer program, to obtain an insight into the general structure of the func- 
tions u p . This will help with a general analytic proof of the subconjecture and 
also provide indispensable information for solving the regular finite initial value 
problem numerically. 

4.3.4 Problems and prospects 

Condition (|56) suggests that data for solutions satisfying the requirements of 
definition (3.1) have to satisfy at i conditions at all orders. Therefore, a gener- 
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alization of the regularity condition to the case of data with non-trivial extrinsic 
curvature requires a general and very detailed investigation of the solutions to 
the constraints near space-like infinity. Such a study has recently been carried 
out (cf. S. Dain, this volume) and the generalization of the analysis of p£| to 
a much larger class of data, which also requires new considerations concerning 
the gauge, is possible. 

The proof of our conjecture, or its generalization to the case of non-trivial 
extrinsic curvature, requires to show the existence of solutions to the finite 
regular initial value problem near / whose underlying manifolds contain a piece 
of the form {|r| < 1,0 < p < po}. Moreover, we need to obtain sufficient control 
on the smoothness of the solutions near J^ Q = {|r| = 1,0 < p < po}. The 
specification of the precise values of k* and q*(k), possibly referring to Sobolev 
spaces, should be part of this proof. In principle it is not only of interest to 
see whether the conjecture can be justified but we would also like to know how 
the logarithmic terms propagate along , if they do so, and what the role and 
meaning of these terms is. 

It would also be desirable to sort out for more general classes of free data 
the type of non-regularities which arise from solving the constraints and from 
solving the evolution equations. Since free data can be quite 'rough' at space- 
like infinity while still being asymptotically Euclidean, there has to be made a 
reasonable choice. There is no point in striving for a generality which is not 
required and of no use in modeling physical systems. However, already for quite 
natural looking free data the constraints can imply terms of the form p k log- 7 p 
on S (cf. S. Dain, this volume). Clearly it is of interest to see to what extent 
our analysis can be generalized to handle such data even if our conjecture and 
its possible generalizations to more general data bear out. It will be the only 
way to assess the generality and the role of the solutions which admit a smooth 
structure at null infinity in the class of all solutions arising from asymptotically 
Euclidean data. 

If the smoothness of the conformal structure obtained at null infinity turns 
out to be much weaker than expected in systems which appear of 'physical 
importance'^], we will have to reconsider definition (3T). While it may not 
be applicable any longer in the desired generality, the underlying idea is viable 
because the conformal field equations may be used to control also solutions with 
a 'rough' asymptotic behaviour (@, cf. also P. T. Chrusciel, this volume). 
If the solutions for suitable data admit polyhomogeneous expansions in r, we 
can expect to relate the coefficients of this expansion, which can be considered 
as functions on again to the initial data. This should allow us to give some 
interpretation for these terms. 

In all such studies the insight into the equations gained by trying to establish 
the extremely sharp characterization of the fall-off behaviour of gravitational 
fields in terms of asymptotic simplicity will prove important. For instance, the 
general analytic proof of the subconjecture stated above will be of particular 

1 Besides getting control on the mathematical possibilities, finding out which data are of 
'physical importance' may turn out to be the most interesting part of the whole exercise. The 
interpretation of initial data is still a wide open field. 
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interest because it will provide a further understanding of the peculiar features 
of the field equations. The precise structure of the lower order terms in the 
equations, which is essentially determined by the fact that we are dealing with a 
soldered gauge theory (cf. the discussion of the general conformal field equations 
in the context of normal conformal Cartan connections in |35fl), becomes critical 
here. Understanding its role in the context of the data dependent conformal 
singularity at space-like infinity may shed light also on the mechanism of the 
field equations in other singular situations. 

The concept of an isolated system is not part of the general theory but a 
matter of expediency (cf. the discussion in |$6|). The final criterion for its 
definition is the possibility to model physical systems in a way which allows us 
to extract interesting information about physical processes. If the conjecture 
above bears out, the resulting space-times will allow us to perform very detailed 
discussions of all their physical aspects and the proof of the conjecture will 
provide information which reduces many discussions to direct calculations. 

An example for this is given by the calculation of the NP-constants in terms 
of the initial data in Jj2| , As neatly discussed in [50], these provide interesting 
restrictions on the global structure of solutions. Another example is the decision 
on the reducibility of the BMS group to the Poincare group in Q . So far both 
results still rely on the correctness of our conjecture. In fact, many earlier 
studies of physical concepts associated with the asymptotic structure will be 
put on a firm analytical footing. 

If logarithmic singularities at null infinity (or the coefficients of the loga- 
rithmic terms) turn out to be of any physical significance this will not only 
create problems for the concept of asymptotic simplicity but also for the con- 
ventional method of calculating gravitational radiation by solving numerically 
initial boundary value problems for Einstein's field equations. In that approach 
the singularities and what might be lost by ignoring them cannot even be seen. 
While the initial data for the initial boundary value problem can be arranged 
not to be affected by the modification of the data near space-like infinity, with 
our present knowledge there appears to be no way, whatever boundary treat- 
ment in the numerical calculation is used, to control whether one is calculating 
a solution where the influence of such a modification enters through the bound- 
ary into the numerical solution (cf. jfi} for the amount of freedom and control 
available in the analytic treatment of the initial boundary value problem for 
Einstein's field equations). 

On the other hand if the latter method makes any sense at all, it may for 
practical purposes be quite sufficient and of advantage to calculate numerically 
asymptotically simple solutions for data with a Schwarzschild end if the pro- 
cedures in pi , pl| are amenable to a numerical treatment. It is quite simple 
to identify in the evolution of such data hyperboloidal hypersurfaces which in- 
tersect J + in the Schwarzschild part of the evolution and which extend with 
the implied data smoothly as solutions to the conformal constraints across J7 + 
(pi]]). In fact, the corresponding hyperboloidal data for the conformal field 
equations can be calculated by solving a system of ODEs ([Q). This would 
resolve a number of difficult problems occurring in the present numerical codes 
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(cf. J. Frauendiener and S. Husa, this volume). 

Again, while this in the end may turn out sufficient for practical purposes, 
we have no possibility to decide on that before the existence of asymptotically 
simple solutions has been discussed under sufficiently general assumptions. In 
that case quite general complete maximal hyperbolic solutions together with 
their asymptotics may become calculable from standard Cauchy data and the 
interpretational framework associated with the concept of asymptotic simplicity 
may allow us to draw the desired conclusion. 

It appears that presently the conformal field equations offer the only way 
to perform such calculations. In the case of the characteristic method, which 
also allows one to calculate numerically semi-global solutions, there may arise 
problems because null hypersurfaces have the tendency to develop caustics (an 
intrinsic feature which cannot be overcome by smoothing procedures, cf. E7| for 
a discussion of caustics in the context of the characteristic initial value problem 
and G. Galloway, this volume, for an illustration of the ubiquity of conjugate 
points on null geodesies). There are also problems in prescribing appropriate 
initial data. Assuming past time-like infinity to be represented by a regular 
point, one could think of prescribing null data on the set J" U which only 
has a very special and well understood 'caustic' at i~ . However, at present it 
is by no means clear how to prescribe these data such that they develop into 
a solution which admits an asymptotically flat Cauchy hypersurface. Similar 
difficulties occur if one tries to invent other ways of calculating complete space- 
times from characteristic data. 

4.4 Time-like infinity 

It is a remarkable feature of the conformal Einstein equations, which should 
not be taken for granted, that they force the null generators of J + to meet 
in precisely one regular point i + of the conformal extension if the solution ex- 
ists for long enough in a sufficiently regular fashion and the initial data are 
prescribed appropriately (tj30|). In such a situation 'time-like infinity' has a 
clear-cut meaning, it is locally similar to that of Minkowski space, and it can 
'readily' be attained in numerical calculations ([Q, cf. also S. Husa, this vol- 
ume). The calculation of i + provides an opportunity to test numerical codes in 
long range calculations of space-times and to perform studies of the non-linear 
interaction of gravitational waves under well controlled conditions, however, it 
is of course not a purpose in itself. 

To calculate the gravitational radiation emitted by the merger of two black 
holes we will have to consider solution which develop singularities and horizons. 
While it may appear wise in a numerical calculation to avoid, if possible, entering 
the horizon and hitting the singularity, such an attitude does not help in an 
analytical investigation and we do, in fact, consider it necessary to attack the 
situation head on. In fluid mechanics there exists a long tradition of studying 
singularities (shocks) and there have been developed methods to handle them in 
numerical calculations. The singularities observed in general relativity, however, 
are in general much more complicated and severe from a conceptional point of 
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view as well as from the point of view of actual analytical expressions. In spite of 
a considerable effort to come to grips with them analytically and, more recently, 
numerically (cf. D. Garfmkle, this volume, and also [jTj and the literature given 
there, which discusses some remarkable achievements), we are still far from a 
situation where we could attain the desired information in sufficiently general 
situations. In fact, the 'well-known' Schwarzschild-Kruskal or Kerr solutions 
still pose difficulties. 

There exist several ways to discuss 'singular points of space-times' in an 
abstract way (cf. jn] and the references given there), however, there are no 
concepts available which would allow us to perform under general assumptions 
an analysis of singularities in the context of the field equations. While it is a 
central problem in singular situations to obtain detailed statements which are 
independent of any specific gauge, it is at least as difficult to develop gauge 
conditions which allow us to arrive at any interesting statement at all. The 
latter problem certainly does not have an universal solution. The best we can 
hope for is to find manageable conditions which by their general nature and by 
experience with specific solutions may be expected to work for certain classes 
of solutions. 

We have seen above that conformal Gauss systems have many nice local 
properties and the discussion of the regular finite Cauchy problem suggests 
that they also behave nicely near space-like infinity. It can also be shown that 
they can be used in calculations of i + from hyperboloidal Cauchy data which 
are sufficiently close to Minkowskian hyperboloidal data. However, in all these 
situations the gravitational field must be considered as weak and one would like 
to know whether conformal Gauss systems can also exist globally and with good 
asymptotic properties on asymptotically flat space-times which posses regions 
of strong curvature or even singularities and horizons. 

Caustics of congruences of metric geodesies arise because the latter obey 
equations of second order. Congruences of conformal geodesies are governed 
by equations of third order (cf. |64fl ) and are thus even more prone to caustic 
formation, in fact, different conformal geodesies may become tangent to each 
other. However, the same structural property which may be responsible for 
such severe caustics allows us to show that there exist conformal Gauss systems 
on the Schwarzschild-Kruskal solution which smoothly cover the entire space- 
time and extend smoothly and without degeneracy through Sf ± into any smooth 
extension of the space-time through null infinity ([0). 

These systems are characterized by initial data for the conformal geodesies on 
the hypersurface S of the Schwarzschild-Kruskal solution which extends analyt- 
ically the hypersurface {t — 0} of the Schwarzschild space. The data considered 
in [ p4| are invariant under the space reflection symmetry in S and the rotational 
symmetry and the congruence is orthogonal to S. On S are given the standard 
spherical coordinates tp, 6 and the isotropic radial coordinate r, < r < oo, 
which is related to the radial coordinate f of the standard Schwarzschild line 
element g s = (1 - 2JE) dt 2 - (1 - ^f-y 1 df 2 - f 2 do 2 on the subset f < 2 m 
by r = l/2(f — m + yjf(r — 2 m)). The two asymptotically flat ends of S are 
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at r = and r = oo. The coordinates <fi, 9, r are dragged along with the con- 
gruence and define together with the parameter r of the conformal geodesies a 
smooth coordinate system (ignoring here the singularity inherent in the spher- 
ical coordinates which in the present situation can be easily avoided) on the 
Schwarzschild-Kruskal space with S = {t = 0}. The initial data are prescribed 
such that the conformal factor takes the form 9 = (r- ! j) 2 (r + ?r)~ 2 (/(?-) 2 -t 2 ) 

where /(r) = r/(r 2 - 

The coordinates stay regular, because of the specific choice of the 1-form b on 
S. The acceleration of the conformal geodesies, which is determined by b, tends 
by our choice to spread out the curves such that no caustics can develop. This 
latter fact is established by analyzing an equation satisfied by the connecting 
vectors of the congruence of conformal geodesies which generalizes the Jacobi 
equation known for metric geodesies. We consider the fact, that such a geometric 
analysis is possible as particularly useful. 

It is well known that the Schwarzschild-Kruskal solution admits global dou- 
ble null coordinates (cf. B. Schmidt, this volume), however, these coordinates, 
or coordinates derived from them by simple algebraic relations, do not extend 
smoothly through ^ and they are very specific for the Kerr family. The prop- 
erties of the conformal Gauss systems mentioned above raise the hope that such 
systems can be constructed in much more general situations. Of course, while 
there are a few general considerations which may help here, a good choice of 
the initial data for the conformal geodesies will depend on the specific solution 
and one has to gather experience with this, as with any choice of gauge. If 
conformal Gauss systems are used in the context of the general conformal field 
equations to perform numerical calculations, the tensor fields entering the gen- 
eralized Jacobi equations are essentially calculated as one goes along and the 
possible tendency of the underlying congruence of conformal geodesies to form 
caustics can be monitored during the calculation. 

It should be noted that Gauss systems for which the initial tangent vector 
field on S and the pull back of the 1-form field b to S are identical, define the 
same congruence of curves, if these are considered as point sets. The parame- 
terizations of these curves, however, may be different and the systems will thus 
define different foliations. This freedom has hardly been explored yet. 

Since the space reflection in S is given by r — > m 2 /4r, it is sufficient to 
discuss the conformal geodesies starting at the subset of S where r > to/2, 
which cover a full f > 2 to part and two 'halves' (one in the past and one in the 
future of S) of the f < 2 m part of the Schwarzschild solution. The conformal 
geodesies with cf) = const., 9 = const., r = r + = (3 + ^/h)m/A, stay on a 
hypersurface {f — 5 to/2} of the Schwarzschild-Kruskal space and approach 
time-like infinity, denoted by i + , as r — > r% = 2/V5to. The conformal geodesies 
starting with r > r + run out to null infinity and arrive there as r — > T scr i = f{r), 
which is consistent with the expression above of the conformal factor. 

The curves starting with ^ < r < r + hit the singularity at f = as 

r -> r sing = f(r)H(K(k),k), where k(r) = V2 (l + ^7j2 (r - f)) ^ 
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K(k) = f£ /2 (l - k 2 sin 2 x)' 1 ' 2 dx, and 



2-k 2 r x dn 2 {v,k)-{l-k 2 ) 
q dn 2 {v,k) + (1 - A: 2 ) 



with Jacobi's elliptic function (in. They pierce the horizon H. + = {r — 2 m} 
at r = r hor = f(r)H(x h ,k), with a* = j o ™(fc 2 /(4-3fc 2 )) (1 _ fe2 ^2^-1/2 dx 
As to be expected, we find that r scr i — > Ti as r — > r + with r > r + and also 
7"sin g , T hor — > n as r — > r + with r < r+. 

Signals periodic in proper time of an observer approaching the horizon are 
received on J + at Bondi-times which grow only logarithmically^. Consequently, 
a high resolution is required in a numerical calculation which aspires to probe 
the field near i + by following outgoing null geodesies which start on a Cauchy 
hypersurface and arrive at null infinity at very late Bondi-times. Therefore we 
consider it particularly important to acquire a detailed understanding of the 
characteristic features of the field near i + . 

It is customary, and natural from the point of view of the causal structure (cf. 
[ffl|), to denote time- like infinity in the schematic (conformally compactified) 
causal diagrams by one or several ideal points i + , but the specific structure of 
the field in the region close to them is little understood. In our representation 
of the Schwarzschild-Kruskal space the behaviour of the term (O f) 2 da 2 of the 
metric Q 2 gs in the limit as t -> r t on the hypersurface {f = f} or {f = 
2 m} suggests that we should consider i + as a point with coordinates r = n , 
r = r + (the coordinates <f> and 9 loosing their meaning there as at the origin 
of a polar coordinate system). Without further analysis, however, this can 
only be considered as a way of speaking. Contrary to what is observed in 
Minkowski space, i + is not only approached by time-like geodesies and by the 
null generators of J + , but also by null geodesies originating from the interior 
of the space-time, namely those generating the event horizon and those which 
stay on the set {f = 3 m} or approach it in the future. Such phenomena do not 
occur for regular conformal structures and thus i + cannot be considered as a 
point in a suitably defined regular conformal extension. 

It is reasonable to expect that conformal geodesies of a regular conformal 
structure which approach a point p in the underlying manifold M can be ex- 
tended, perhaps after an admissible reparametrisation, through that point as a 
smooth solution to the conformal geodesic equation, and that the curve will be 
time-like at p if it was time-like before (while such a statement is known to hold 
for metric geodesies it still has to be shown for conformal geodesies) . Assuming 
this to be true, we can observe an effect of the conformal singularity near i + 
in our global conformal Gauss system. Since dr scr i/dr converges to a finite 
negative number as r — ► r + with r > r + , we see that J + approaches, as to be 
expected, the hypersurface {r = r + } near i + under a certain angle. In contrast, 

vg-l 

we find that dr S i ng /dr fa (r + — r) ^ as r — * r + with r < r+, such that the 

2 I owe much to discussions with P. Hiibner about this point. However, we both have to 
acknowledge the priority of W. Allen's formula: eternity is very long, especially near the end. 
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graph of T S i n g, and with it also the graph of Th or , becomes tangent to {r = r + } 
at i + . As a consequence, the conformal geodesies generating {r = r + } become 
tangent to the null geodesies generating H + when they approach i + . 

Though this may look like a subtlety, one may wonder about its consequences 
if one tries to observe the CFL condition in a numerical calculation near i + . 
Moreover, the speed with which the singularity and the horizon approach the 
hypersurface {r — r + } may require special measures to sufficiently resolve these 
structures numerically. 

There are of course much more drastic effects. While the rescaled confor- 
mal Weyl tensor (expressed in suitably adapted coordinates or frame fields) 
extends smoothly to J + , it diverges along the null generators of the horizon 
like _1 near i + , and blows up like 1/f 3 near the singularity (if the conformal 
factor approaches a positive value there, as it does in our case), resulting in a 
non- uniformity with strong gradients near i + . Moreover, slight changes in the 
data may lead to a Kerr solution with drastically different behaviour near i + , 
including a Cauchy horizon, a completely different singularity structure, and 
causality violations (cf. Q). In recent years considerable work has been done 
to understand what happens to these structures under perturbations of the data 
(cf. Q and the references given there) but these studies supply surprisingly 
little precise information on the behaviour of the field in the close vicinity of i + . 

One may object to our discussion that we create unnecessary problems by 
insisting on conformal rescalings in a region where the field is usually considered 
as weak. But the phenomena we discuss are present anyway, they are just 
brought out most succinctly by suitable rescalings. ft is not so clear, though, 
what 'suitable' may mean here. In the case of space-like infinity we have seen 
that there are representations which allow us to see properties which are hard 
to analyse with other choices of the gauge. 

This raises the question whether there exists a kind of 'resolution of the 
singularity at time-like infinity' which makes the situation more amenable to a 
detailed analysis and to a numerical treatment. The situation near i + is clearly 
quite different from that near i°. There is nothing near i which compares 
with the singularity '{f = 0}' and the horizon TL + and the way they meet with 
J + at i + . Because of these structures the singularity obtained in a conformal 
representation which makes time-like infinity a point is much less isotropic than 
the singularity obtained if space- like infinity is represented by a point. The 
singularity at f = cannot be removed by a simple rescaling without shifting 
it to infinity. Moreover, i + is causally differently related to the physical space- 
time than i°. A 'resolution of the singularity at i + ' must be compatible with 
a time evolution which essentially takes into account the complete information 
given by the initial data. It should be emphasized that we are not interested 
here in a formal definition of a space-time boundary which relies on more or less 
implicit assumptions of the time evolution of the field. We seek a construction 
which will help us analyse the latter. This construction can only be fruitful if it 
follows from the structure of the field equations in a similar way as the cylinder 
at space-like infinity does. 

Looking again at the complicated behaviour of the fields near time-like infin- 
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ity suggested in Q, this task may appear hopeless. However, it seems to be the 
accepted view that in the case of a gravitational collapse to a black hole some 
kind of universal behaviour can be expected near i + in the sense that the field 
will in general approach near time-like infinity somehow that of a Kerr solution 
(p3)). The recent proof of the Penrose inequality (pS, Q) lends some support 
to this, because the original argument leading to this inequality relied on the 
'establishment viewpoint' indicated above (p9fl). 
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